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, The band structure of the Lame equation, viewed as a one-dimensional Schrodinger 

equation with a periodic potential, is studied. At integer values of the degree pa- 
rameter £, the dispersion relation is reduced to the £ — 1 dispersion relation, and a 
previously published £ — 2 dispersion relation is shown to be partly incorrect. The 
Hermite-Krichever Ansatz, which expresses Lame equation solutions in terms of 
£ — 1 solutions, is the chief tool. It is based on a projection from a genus-^ hyper- 
I elliptic curve, which parametrizes solutions, to an elliptic curve. A general formula 

J> ■ for this covering is derived, and is used to reduce certain hyperelliptic integrals 

lO I to elliptic ones. Degeneracies between band edges, which can occur if the Lame 

■ equation parameters take complex values, are investigated. If the Lame equation is 
viewed as a differential equation on an elliptic curve, a formula is conjectured for 

■ the number of points in elliptic moduli space (elliptic curve parameter space) at 
I which degeneracies occur. Tables of spectral polynomials and Lame polynomials, 
i i.e., band edge solutions, are given. A table in the older literature is corrected. 
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1. Introduction 

(a) Background 



' The term 'Lame equation' refers to any of several closely related second-order or- 

I dinary differential equations ijWhittaker fc Watso'nlll927t lErdelvilll953-55t lArscottI 

11964) . The first such equation was obtained by Lame, by applying the method 
of separation of variables to Laplace's equation in an ellipsoidal coordinate sys- 
tem. Lame equations arise elsewhere in theoretical physics. Recent application ar- 
eas include (i) the analysi s of preh eating af ter inflation, arising from paramet- 
ric am plification l(Bo vanov skv et alJll9 96: G reene et al.lll99 7: Kaiser 1998; Ivano^ 
ii) the stabil ity analysis of critical droplets in bounded spatial domains 
ijMaier fc Steir 1200 ll): (iii) the stab ility analysis of static configurations in Joseph- 



son junctions i Caputo et al.ll200ol) . (iv ) the computation of the distance-redshift 



relation in inhomogeneous cosmologies llKantowski fc T homas 2001), and (v) mag- 
netostatic problems in triaxial ellipsoids ()Dobne'nr^itteri.l998i) . 

In some versions of the Lame equation, elliptic functions appear explicitly, and 
in others (the algebraic versions) they appear implicitly. The version that appears 
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most often in the physics hterature is the Jacobi one, 

d2 



, „+£(£+ l)msn^(a\m) 



* = (1.1) 



which is a one-dimensional Schrodinger equation with a doubly periodic potential, 
parametrized by m and £. Here sn(-|TO) is the Jacobi elliptic function with modular 
parameter m. m is often restricted to (0, 1), though in general m S C \ {0, 1} is 
allowed. When m G (0, 1), the function sn^(-|m) has real period 2K :— 2K(m) and 
imaginary period 2iK' := 2iK(l — to), with K(to) the first complete elliptic integral. 

If a is restricted to the real axis and m and £ are real, Hl.l() becomes a real- 
domain Schrodinger equation with a periodic potential, i.e., a Hill's equation. Stan- 

dard results on Hill's equation apply ( Magnus & Winkler 1979; McK ean & van Moerbekj 
Il979|) . Equipping 1)1.1(1 with a quasi-periodic boundary condition 

*(a + 2K) ^ e"=(2^%(a) :^ ^ *(a), (1.2) 

where fc G M is fixed, defines a self-adjoint boundary value problem. For any real k 
(i.e., for any Floquet multiplier ^ with |^| = 1), there will be an infinite discrete set 
of energies S R for which this problem has a solution, called a Bloch solution with 
crystal momentum k. Each such E will lie in one of the allowed zones, which are 
intervals delimited by energies corresponding to ^ = ±1, i.e., to periodic and anti- 
periodic Bloch solutions. These form a sequence Eq < Ei < E2 < < E4 < ■ ■ ■ , 
where Eq is a 'periodic' eigenvalue, followed by alternating pairs of anti-periodic 
and periodic eigenvalues (each pair may be coincident). The allowed zones are the 
intervals [E2j, E2j+i]. The complementary intervals (i?2j+i, i?2j-(-2) are forbidden 
zones, or lacuna. Any solution of Hill's equation with energy in a lacuna is unstable: 
its multiplier ^ will not satisfy |£| = 1 , and its crystal momentum k will not be real. 

It is a celebrated result of lined that if the degree £ is an integer, 

which without loss of generality may be chosen to be non-negative, the Lame 
equation (|l.l|l will have only a finite number of nonempty lacunae. A converse 
to this statement holds as well i Gesztesv fc Weikardlll995(7l) . If £ is an integer. 



the Bloch spectrum consists of the £+1 bands [£^0, ^'i], [E2, E3], . . . , [E2£, 00), and 
£{£ + l)TOsn(-|TO) is said to be a finite-band or algebro-geometric potential. The 
2£+l band edges Eq, . . . , E21 are algebraic functions of the parameter m. That is, 
they are the roots of a certain polynomial, the coefficients of which are polynomial 
in TO. The corresponding periodic and anti-periodic Bloch solutions are called Lame 
polynomials: they are polynomials in the Jacobi elliptic functions sn(a|m), cn(Q!|TO), 
and dn(a|m). The double eigenvalues embedded in the topmost band [i?2£,oo) (the 
'conduction' band), namely E2j = £^2^+1, j > £, are loosely called transcendental 
eigenvalue s. For £ = 1 at leas t, they are known to be transcendental functions of to 
ijChud nov ksv fc Chu dnovsk^ ll98(i . 

There has been much work on algebraizing the integer-^ Lame equation, to 
facilitate the computation of the band edges and the coefficients of the Lame poly- 
nomials (Alhassid et al. 1983; Turbiner 1989; Li & Kusnezov 1999; Li et al. 200^ 
iFinkel et alJl200Clfl . Such schemes have been extended to the case when £ is a half- 
odd-integer, in which there are an infinite number of lacunae. In this case, certain 
'mid-band' Bloch functions, namely ones with ^ = ±i and real period 8K, are al- 
gebraic functions of sn(Q;|TO). Certain rational values of £ with 2£ ^ 7L also yield 
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algebraic Bloch functions, provided the parameters m and E are chosen appropri- 
ately ilvlaicr 2004). 

An algebraic understanding of band edges is useful, but it is also desirable to 
have a closed-form expression for the dispersion relation: fc as a function of E. The 
value of k is not unique, since it can be negated (equivalently, ^ i-^ ^^'^ 
integer multiple ofn/K can be added. However, each branch has the property that 
k ^ E^l'^ or ~ —E^l"^ to leading order s& E ^ -|-oo. Also, fc € R in each band. 

The goal of this paper is the efficient computation of the dispersion relation 
when t is an integer. The following are illustrations of why this is of importance 
in theoretical physics. In application (i) above (preheating after inflation), particle 
production is due to parametric amplification: a solution having a multiplier ^ with 
1^1 > 1. This corresponds to the energy E not being at a band edge, or even in a 
band, but in a lacuna. In application (ii) (the stability analysis of a critical droplet), 
the analysis includes an imposition of Dirichlc t rather than qu asi-periodic boundary 
conditions on an f = 2 Lame equation ifMaier fc SteinllioOlf) . The resulting Bloch 
solution is not a Lame polynomial, but rather a mid-band solution. 

When i is an integer, the Lame equation is integrable, and the general integral 
of can be expressed in terms of Jacobi theta functions. The dispersion relations 
k — ki(E\m), € > 1, can in principle be computed in terms of elliptic integrals. The 
case £ = 1 is by far the easiest. If ^ = 1, the solution space of except when 

i? is at a band edge, will be spanned by the pair of Hcrmite-Halphen solutions 

$(a; ±ao|m) ^ T^T^ Z(±ao|m)] . (1.3) 

Here H, 8, Z are the Jacobi eta, theta and zeta functions, with periods 4i^, IK, 2K 
respectively, and ao € C is defined up to sign by dn^(Q!o|™) ^ E — m. So 

ki{E\m) ^ -iZ{ao\m) +T:/2K{m), dn^{ao\m) = E - m, (1.4) 

up to multi-valuedness. This is a parametric dispersion relation. It has been ex- 

ploited in a study of W annier-Stark resonances with non-real E and k ()Grecchi fc Sacchettil 
Il997t ISacchettilll997^ . However, the extension to ^ > 1 is numerically nontrivial. 
k£{E\m) turns out to equal J2r=o[~^'^i^r\iTi) -I- 7r/2K (TO)], where |arlf~n satisfy 
coupled transcendental equations involving E and m l|Whittaker fc Watsonlll927[ 
§23.71). Li, Kusnezov fc lachello calculated and graphed fc = k^jElm) as well a s 
k = ki{E\m) in the 'lemniscatic' case m = ^ l|Li fc Kusnezovlfl999l : iLi et al.ll2000(l . 
Unfortunately their graph of fc = k2{E\^) is incorrect, as will be shown. 

(b) Overview of results 

When £ > 1, we abandon the traditional Hermite-Halphen solutions, and exam- 
ine instead the implications for the Lame dispersion relations of what is now called 
the Hermite-Krichever Ansatz. This is an alternative way of generating closed-form 
solutions of the Lame solution at arbitrary energy E; for small integer values of £, 
at lea st. Until the 1980s, the only reference for the Ansatz was the classic work of 
iHalphe n ( 1888, c hapitr e XII), who applied it to the cases £ = 2,3, 4, and in part to 
£ = 5.|Krichev^ l)l98Cl|) revived it as an aid in the construction of elliptic solutions 
of the Kortewe g-de Vries and other inte grable evolution equations, feclokolos et all 
l(l98fil) and Bel okolos fc Enol'skiil 1 200(|) summarize early and recent developments. 
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The Hermite-Krichever Ansatz is easy to explain, even in the context of the 
Jacobi form of the Lame equation, which is not the most convenient for symboHc 
manipulation. It asserts that for any integer £ > 1, fundamental solutions of 
can be constructed from the £ — 1 solutions ±Q!o|rn) as finite series of the form 



exp(±KfQ;), (1.5) 



where the parameter ao is now computed from a reduced energy Ei by the formula 
dn^(ao|'7^) = £e — m. The reduced energy Sg, the exponent k^, and the coefficients 
{cj^^} will depend on E and m. ££{E\m) may be chosen to be rational in E and m, 

and ±K£ to be of the form k£{E\m) times ±i-y/ L({E\m), where ki{E\m) is also 

rational, and Lg{E\m) is the spectral polynomial Y[lLo[-^~ a degree-(2^+l) 
polynomial in E the coefficients of which, as noted, are rational in to. 

If the Lame equation can be integrated in the framework of the Hermite- 
Krichever Ansatz, it follows from H1.5|l that up to sign, etc., 

ke{E\m) = ki{£e{E\m)\m) + ki{E\m)^jLiiE\m). (1.6) 

The dispersion relation for any integer £ can be expressed in terms of the 1 = 1 
relation. To compute fc^, only one transcendental function (i.e., ki) needs to be 
evaluated, since the other functions in H1.6|l are elementary. The only difficult matter 
is choosing the relative sign of the two terms, since each is defined only up to sign. 

The functions £i{E\m), ki{E\m) are rational with integer coefficients, but work- 
ing them out when £ is large is a lengthy task. In principle one can write down a 
recurrence relation for the coefficients {Cj}, and work out £g{E\m), kt{E\m) from 
the condition that the series terminate. However, their numerator and denominator 
degrees grow quadratically as £ increases. This explains why Halphcn's treatment 
of the f = 5 case was only partial. In a series of papers, Kostov, Enol'skii, and col- 



laborators used computer 


algebra systems to perform a full analysis of the cases £ = 


2. 3. 4. 5 JGerdt & Kostov 1989 


Kostov & Enol'skii 1993; 


Enol'skii & Kostov 


1994 



19941) . When £ = b, using Mathematica to compute the inte- 



ger coefficients of rational functions equivalent to £i{E\m), ki{E\m) required seven 
hours of time on a Sparc-1, a Unix workstation of that era {Eilbcck & Enol'sldl 
Il994l) . Until now, their analysis has not been extended to higher £. 

While performing extensive symbolic computations, we recently made a discov- 
ery which is formalized in Theorem L below. For all integer £ > 2, the degree-£ Lame 
equation can he integrated in the framework of the Hermite-Krichever Ansatz, and 
the rational functions that perform the reduction to the £ = 1 case can he com- 
puted by simple formulas from certain spectral polynomials of the degree-£ equa- 
tion, which are relatively easy to work out. These are the ordinary spectral polyno- 
mial Hslot^ ~ Eisim)] associated with the band-edge solutions, and the spectral 
polynomials associated with two other types of closed-form solution that have not 
previously been studied in the literature. We call them twisted and theta-twisted 
Lame polynomials. In the context of the Jacobi form, the former are polynomials in 
sn(a|m), cn(a|TO) and dn(a|TO), multiplied by a factor exp(Ka). (If k € R, 'canted' 
would be better than 'twisted'.) The latter contain a factor resembling (|1.3|l . 
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Theorem L follows from modern finite-band integration theory: specifically, from 
the Baker- Akhiezer uniformization of the relation between the energy and the crys- 
tal momentum. This uniformization is closely tied to classical work on the Lame 
equation (the parametrized Baker-Akhiezer solutions of the integer-^ Lame equa- 
tion are in fact equivalent to the Hermite-Halphen solutions). Theorem L greatly 
simplifies the computation of higher-^ dispersion relations. It also has implications 
for the theory of hypere lliptic reduc tion : the reduction of hyperelliptic integrals to 
elliptic ones JBelokolos et al.lll98(i) . This is on account of the following. For any 
£ > 1 and m, the solutions of the Lame equation, both the Hermite-Halphen so- 
lutions and those derived from the Hermite-Krichever Ansatz, are single- valued 
functions not of E, but rather of a point {E, v) on the -^'th Lame spectral curve: 
a hyperelliptic curve comprising all {E, D) satisfying 

i^' = \{[E-EM)]- (1-7) 

The fact that v = v{E) is two- valued (except at a band edge) is responsible for the 
two-valuedness of, e.g., the parameter — ao(i?) of (|1.3|l . and in general, for the 
uncertainty in the sign of k. The £'th spectral curve generically has genus £ and may 
be denoted Vi := Viijn). For any integer t > 2, the E i— > £i{E\m) reduction map 
of the Hermite-Krichever Ansatz induces a covering tt^ : — > Fi. The first known 
covering of an elliptic curve by a higher^enus hyperelliptic curve was constructed 
by Legendre and generalized by Jacobi {Bclokolos ct al. 1986, § 2). But it is difficult 
to enumerate such coverings, or even work out explicit examples. Those generated 
by the Ansatz applied to the Lame equation are a welcome exception. 

The integral with respect to E of any rational function of E and •\/n(i?), where 
n is a polynomial, is a line integral on the algebraic curve defined by = n(ii^). 
The covering tt^ : F^ — > Fi, a formula for which is provided by Theorem L, reduces 
certain such hyperelliptic integrals to elliptic ones. In modern language, the theorem 
specifies how certain holomorphic differentials on hyperelliptic curves can arise as 
puUbacks of holomorphic differentials on elliptic curves. 

We investigate the degeneracies of the band edges {i?s(m)}^^Q that can occur 
when th e modular pa rameter m is non-real. Such level-crossings were first consid- 
ered by ICohnI l)l888|) in a dissertati on that se ems not to have been followed up, 
though it was later cited bvfwhitta ker fc Wat son (1927, §23.41). We conjecture a 
formula for the ^-dependence of the number of values of rn G C \ {0, 1}, or equiva- 
lently the number of values of the Klein invariant J e C, at which two band edges 
coincide. When this occurs, the genus of the hyperelliptic spectral curve F^ is re- 
duced from £ to — 1, though its arithmetic genus r emains equal to I. Band-edge 
degeneracies are responsible for a fact discovered bv iTurbined l)l989j) : if ^ > 2, the 
complex curve comprising all points {(m, £^s(to))}s1q in C \ {0, 1} x C has only 
four, rather than 2^ -t- 1, connected components. 

This paper is organized as follows. Section |21 introduces the Lame equation 
in its elliptic-curve form and relates the Hermite-Halphen solutions to the Baker- 
Akhiezer function. In §|21 Lame polynomials in the context of the elliptic-curve form 
are classified. In §01 the Hermite-Krichever Ansatz is introduced, and Theorem L 
is stated and proved. The application to hyperelliptic reduction is covered in §|21 In 
§§iniand|3 dispersion relations are worked out and the abovementioned dispersion 
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relation for the case ^ = 2 is corrected. The £ = 3 dispersion relation is graphed 
as well. Finally, an area for future investigation is mentioned in §|S1 



2. The elliptic-curve algebraic form 

In §§|31through|Slwe use exclusively what we call the elliptic-curve algebraic form of 
the Lame equation, which is the most convenient for symbolic computation. In this 
section we derive it, and also define a fundamental multi-valued function $, which 
appears in the elliptic-curve version of both the Hermite-Halphen solutions and the 
Hermite-Krichever Ansatz. 



(a) An elliptic- curve Schrddinger equation 

Many algebraic forms can be obtained from l|l.l|l by changing to new inde- 
pendent variables which a re elliptic function s of a , such as sn(Q;|m). (See, e.g., 
lArscott k, Khabazal l)l962l § 1.1) and lArscottI l)l964 pp 192-3)). A form in which 
the domain of definition is explicitly a cubic algebraic curve of genus 1, i.e., a 
cubic elliptic curve, can be obtained as follows. First, the Lame equation is re- 
stated in terms of the Weierstrassian function p — p{u; g2, gs)- This is the canon- 
ical elliptic function with a double pole at m = 0, satisfying (p')^ = /(p) where 
f{x) := Ax^ — g2X — gs = ^Yl^^iix — e^). For elhpticity the roots {e^ji^^^ must 
be distinct, which is equivalent to the condition that the modular discriminant 
A := g2 — 27(/| be non-zero. Either of 52, 53 G C may equal zero, but not both. 

The relation between the Jacobi and Weierstrass elliptic functions is well-known 
l Abramowitz fc Stegunlll965[ §18.9). Choose {e^}'l^i according to 



,9/'2 — m 2rn — 1 —(m + l)\ 
(ei,e2,e3) = AM ^— ^ , (2.1) 



where A G C \ {0} is any convenient proportionality constant. Then 

_^ 4(m^-m+l) _ g 4(m-2)(2m-l)(m-|-l) 

^ ' y3 2y ' v ■ / 

and the dimensionless (A-independent) Klein invariant J := .gf/^ will be given by 



4 (to^ - m + 1)^ 
27 m^(l — m) 

The two sorts of elliptic function will be related by, e.g.. 



J-^ 'l.n 1^ ■ (2-3) 



an\Az\m) = ns''{Az\m) = ^ ^, (2.4) 

p(z) - 63 ei - 63 

and the periods of p, denoted 2uj, 2lo' , will be related to those of sn^ by 

2uj = 2KIA, 2uj' = 2iK'IA. (2.5) 

The case when 2K, 2K' are real, or equi yalentlv a; € R. o;^ € iR (we a ssume A e R) , 
is the case when 52,53 € R and A > l Abramowitz fc Stegunlll965[ § 18.1). 
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Choosing for simplicity ^ = 1, so that ei — 63 = = 1, and rewriting the 
Lame equation with the aid of H2.4|l . yields the Weierstrassian form 

- + 1)p(m;52,53) + -5 = 0, (2.6) 

where + \K' . (The translation of by \K' replaces msn^ by ns^.) Here 

B:=-^(ei-e3)-f(£+l)e3, i.e., 

B :=-£;+ i^(£+ l)(m+ 1), (2.7) 

is a transformed energy parameter. Changing to the new independent variable x := 
p(u; (72,53) converts (12.61) to the commonly encountered algebraic form 



] dx^^ 2^ X- 
\ 7=1 



d £{£+l)x + B 
4117=1(^-67) 



* = 0. (2.8) 



This is a differential equation on the Riemann sphere := C U {00} with regular 
singular points at a; = ei, 62, 63, 00. Any solution of the original Lame equation 
or the Weierstrassian form H2.6|l which is quasi-periodic in the sense that it is 
multiplied by ^, ^' G C \ {0} when a ^ a + 2K or a ^ a + 2\K' respectively 
(equivalently, when u <— u -\- 2lj or u ^ m + 2cj'), will be a path-multiplicative 
function of x. That is, it will be multiplied by ^, ^' when it is analytically continued 
around a cut joining the pair of points a; = 62, 63 or x = ei, 62, respectively. In the 
context of the algebraic form, the dispersion relation is still a relation between 
the energy and a multiplier ^, but the multiplier is interpreted as specifying not 
quasi-periodicity on C, but rather multi-valuedness on P^. 

The algebraic form (|2.8() of the Lame equation lifts naturally to the complex 



elliptic curve Eg^.g^ := {{x,y) G | = f{x)} U {(00,00)} over P^, parametrized 
by (cc, y) := {p{u), p'{u)). One may rewrite H2.8(l in the elliptic-curve algebraic form 

[e{i+l)x + B]^^' ^0. (2.9) 

This is an elliptic-curve Schrodinger equation of the form 

, \ 2 

' ^ 5- = -B^-, (2.10) 




^dx) '^'^^^^ 

with the (rational) potential function q{x) taken to equal £{£ -\- l)x. Equation (|2.9|l 
follows directly from (|2.6|) . since d/du = p'd/dp = yd/dx. It is a differential 
equation on Eg^^g^ with a single singular point: a regular one at {x,y) = (00,00). 
Note that the 2-to-l covering map tt : Eg^^g^ P^ defined by TT{x,y) = x has 
{(e^, 0)}i^^]^ and (00,00) as simple critical points. One reason why (|2.9|l is more 
fundamental than (|2.8|l is that the singular points of l|2.8|) at x = ei, 62, 63 can be 
regarded as artifacts: consequences of {(e^, 0)}i^^;^ being critical points of tt. 

The complex-analytic differential geometry of the elliptic curve Eg^ ^^ takes a bit 
of getting used to. Both x and y are meromorphic P^-valued functions on Eg^ ^^, 
and the only pole that either has on Eg^_g^ is at the point O := (00,00). In a 
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neighborhood of any generic point (a;, y) other than O and the three points (e-y, 0), 
either x ot y will serve as a local coordinate. However, near each (e^, 0) only y will 
be a good local coordinate, since dy/dx diverges at x = e^. Also, x has a double 
and y has a triple pole at O, so the appropriate local coordinate near O is the 
quotient x/y. The 1-form dx/y is not merely meromorphic but holomorphic, with 
no poles on Eg^^g^. Its dual is the vector field (or directional derivative) yd/dx. 

Elliptic functions, i.e., doubly periodic functions, of the original variable u £ C 
correspond to single- valued functions on Eg^^g^. These are rational functions of x,y 
and may be written as Ro{x) + Ri{x)y, i.e., Ro{p{u)) + Ri{p{u))p' (u) . The formula 
{yd/dx)y = 6a;^ — ^52 allows such functions to be differentiated algebraically. In a 
similar way, quasi-doubly periodic functions of u (sometimes called elliptic functions 
of the second kind), which are multiplied by ^ when m <— u + 2w and by ^' when 
u <— u + 2lu' , correspond to multiplicatively multi- valued functions on Eg^^g^. 

^92 93 genus 1 and is topologically a torus. A fundamental pair of loops 
that cannot be shrunk to a point may be chosen to be a loop that extends between 
(e2,0) and (63, 0), and one that extends between (ei,0) and (e2,0), with (if (?2,53 
are real, at least) half of each loop passing through positive values of y, and the 
other half through negative values. One way of constructing an elliptic function 
of the second kind is to anti-differentiate a rational function R(x,y). Here 'anti- 
differentiate' means to compute / R{x,y) dx/y, its indefinite integral against the 
holomorphic 1-form dx/y. The resulting function will typically have a non-zero 
modulus of periodicity associated with each loop. If so, exponentiating it will yield 
a multiplicatively multi-valued function on Eg^^g^, with non-unit multipliers 



(b) The Hermit e-Halphen solutions 



The fundamental multi-valued function $ on the elliptic curve will now be de- 
fined. It is an elliptic-curve version of Halphen's I'element simple (Halphcn 1888). 

Definition. On the elliptic curve Eg^^g^ , the multi-valued meromorphic function <I>, 
parametrized by (a;o, yo) G Eg^^g^ \ {(00, 00)}, is defined up to a constant factor by 
a formula containing an indefinite elliptic integral. 



^ix,y;xo,yo) = exp 



- [ [ y^y° \ ^ 

2 J \x- xoj y 



(2.11) 



Its multi-valuedness, which is multiplicative, arises from the path of integration 
winding around Eg^^g^ in any combination of the two directions. Each branch 
of $ has a simple zero at {x, y) — {xq, yo) and a simple pole at [x, y) — {00, 00). 

To motivate the definition of $, a brief sketch will now be given of the construc- 
tion of the Hermite-Halphen solutions of the elliptic-curve algebraic Lame equa- 
tion H2.9|l for integer £ > 1. The standard published exposition is not fully algebraic, 
being framed largely in the context of the Weierstrassian form ijWhittaker fc WatsorJ 
I1927I §23.7). This sketch will relate the Hermite-Halphen solutions to modern 



finite-band integrati on theory and the concept of a Baker- Akhiezer function ijCesztesv fc HolderJ 
l2nn3HTreibichll200ll) . The starting point is the differential equation 



dx 



- 4[q{x) + B](y—]-2 



dx 



q{x) 



■ = 0. 



(2.12) 
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Table 1. Hermite-Halphen polynomials tvan der Waall\200i. table A. 2) 



£ Ti{x;B;g2,g'j.) 



1 x-B 



2 



The differential operator in (|2T2|) is the 'symmetric square' of the elliptic-curve 
Schrodinger operator of H2T()|I , so the solutions of l|2T2|l include the product of any 
pair of solutions of (|2TU|) . If the potential function q{x) is rational, it is known that 
the solution space of (|2T2I) contains a function T{x\ B) which is (i) meromorphic 
in X (the only poles being at the poles of q{x)) and (ii) monic polynomial in _B, 
if and only if H2T0|1 is a finite-band Schrodinger equation on E„^ , i.e.. a finite -band 
Schrodinger equation with a doubly periodic potential ("its & Matveev"l974^. This 
is an alternative to the characterization of Gesztes y & Wcik ard (199(3), according 
to which H2.10|l is finite-band if and only if for all B € C, every solution ^> is 
meromorphic on Eg^^g^ (multi-valuedness being allowed). 

For example, when q(x) := £{£+l)x, there is a polynomial solution J-e{x; _B; 52, 53) 
of H2.12(l that satisfies conditions (i) and (ii), of degree i in both x and B. So the 
integer-^ Lame equation is a finite-band Schrodinger equation. !Fi{x;B]g2,gs,) is 
called the ^'th Hermite-Halphen polynomial. It may be written J^i{x; B; g2, gs) 
when it is normalized to be monic in x, rather than in B (see table 

In the case of a general rational potential q{x) which is finite-band, let J^{x; B) 
denote the specified solution of H2.12|l . and let its degree in B be denoted £. It follows 
from manipulations parallel to those of Whittaker & Watson that the function 
on Eg^,g^ defined by a formula containing an indefinite integral, 



will be a solution of the Schrodinger equation (|2.10() . Here B E C and is a in- 
dependent but position-independent quantity, determined only up to sign, that can 
be computed by what Whittaker & Watson term an "interesting formula," 



(It is not obvious that the right-hand sid e is independent of the point {x,y) G 
^32,33-) It is widely known ijSmirno 3|200i) that ly is identical to the coordinate v 
on the spectral curve defined by — Yll^^QiB — Bs), where {i3s}slo the band 
edges of the Schrodinger operator; though no really simple proof of this fact seems to 
have been published. A consequence of this is that the formula (|2.13l) parametrizes 
solutions of the elliptic- curve Schrodinger equation (12.11)11 bv (B, v) 6 F^\{(cx), cx))}. 
As defined, 4" is called a Baker-Akhiezer function terichevedllQQfil) . 

Consider now the special case of the integer-^ Lame equation H2.9|l . In this 
case the function ^E* computed by (|2.13|l from the Hermite-Halphen polynomials 
J-^ = T(,{x; B; g2, 53) is in fact an Hermite-Halphen solution of the Lame equation, 
reexpressed in terms of the elliptic curve coordinates {x,y). One can write 4" = 
^^{x,y; B; g2, gs), where the superscript '±' refers to the ambiguity in the sign of 
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V = viS). If z/ 7^ 0, the two solutions are distinct. They are path- multiplicative, 
since they are exponentials of anti-derivatives of rational functions on £"^2, ga- 
it should be noted that the Hermite-Halphen polynomials are not merely a 
tool for generating the soluti ons of th e Lame equation. They are algebraically 
interesting in their own rig ht. lKleinI l|l892l figures 1,2) supplies a sketch of the real 
portion of the curve Ti{x\ B) = when ^ = 5,6, showing how when £ > A, B is a 
band edge only if Ti{x; B), regarded as a polynomial in x, has a double root. 

The relevance of the fundamental multi- valued function $ can now be explained. 
It follows from H2.13|l . and the fact that Ti = B — x (see table that 

^afix, V- B- 52, <?3) = y; ^,±^453 - g^B - g^). (2.15) 

That is, if {xqjUq) G Eg^^g^ is 'above' = B, then <&(•, -^xq, j/o) will be a solution of 
the t = 1 Lame equation in the form H2.9|l . There are two such points, related by j/o 
being negated, unless 4i?3 — g^B — ^3 = 0, i.e., unless B = Ci, 62, 63, in which case 
2/0 = is the only possibility. These are the three band-edge values of B for 1=1. 

It is not difficult to show that the £ = 1 Hermite-Halphen solutions H2.15|l 
are identical to the solutions H1.3|) . though they are expressed as functions of the 
variable {x^y) £ Eg^^g^ rather than the original independent variable a € C. The 
parametrizing point (xo,±yo) G Eg^^g^ corresponds to the parameter ±ao € C 
of H1.3|l . These solutions are clearly easier to formulate in the elliptic curve context. 

For any integer the Lame dispersion relation can be computed numerically 
from (|2.13|) by calculating the multiplier arising from the path of integration wind- 
ing around £"92,33- However, (|2.13l) is not adapted to symbolic computation. By 
expanding the integrand in partial fractions one can derive the remarkable formula 

e 

■i>f{x,y;B;g2,g3) = JJ^x^y; Xr,yt), (2.16) 

r=l 

where {{xr,y^)}i^i are points on Eg^^g^ abo ve ^Xr}i-^ , the -B-dependent roots of 
the degree-^ polynomial Ti{x; B; g2, gs). (Ci. IWhittaker fc WatsonI l)l927l §23.7).) 
Unfortunately, when £ > 5 the roots {xr}r=i cannot be computed in terms of 
radicals. This reduction to degree-1 solutions is less computationally tractable than 
the one that will be provided by the Hermite-Krichever Ansatz. 



3. Finite families of Lame equation solutions 

The solutions of the integer-^ Lame equation include the Lame polynomials, which 
are the traditional band-edge solutions. In the Jacobi-form context they are periodic 
or anti-periodic functions on [0,2/^], with Floquet multiplier ^ = ±1, respectively. 
There are exactly 2£+l values of the spectral parameter B G C, i.e., of the energy E, 
for which a Lame polynomial may be constructed, the counting being up to multi- 
plicity. By definition these are the roots of the spectral polynomial Lg{B; g2, gs). 

As functions on the curve Eg^ ^g^ , the Lame polynomials are single or double- 
valued and are essentially polynomials in the coordinates x, y. (In the Weierstrassian 

context p, p' substitute for x,y.) However, no fully satisfactory table of the Larne 

polyn omials or the Lame spect ral polynomials h as yet been published. Whittaker fc Watso: 
l(l927i § 23.42) refer to a fist of lGiierritor'3 l(l909>) that covers £ < 10. Sadly, although 
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he produced it as a dissertazione di laurea at the Univers ity of Naples , most of his 
results on ^ > 5 are incorrect. This has long be en known llStruttlll967l), b u t his pa- 
per is still occasionally cited for completeness (|Gesztesv k, Holdenll2003|) . lArscottI 
1)19641 § 9.3.2) gives a brief table of the Jacobi-form Lame polynomial s, covering 
only l = 1,2,3. His tab le is correct, with a single misprint jFernand ez C. et al.l 
l200Ct iFinkel et al.ll200ol) . But its brevity has been misinterpreted. An erroneous 
belief has arisen that when ^ > 4, the Lame polynomial coefficients and band edge 
energies cannot be expressed in terms of radicals. This sets in only when ^ > 8. 

Due to these confusions, in this section we tabulate the Lame polynomials and 
the spectral polynomials L^(i3; 172, ffs). Both are computed from coefficient recur- 
rence relations. We supply such relations and tables of spectral polynomials for 
the twisted and theta-twisted Lame polynomials, as well. The number of values 
of -B G C for which the latter two sorts of solution exist, i.e., the degrees of their 
spectral polynomials, will be given. All three sorts of solution will play a role in The- 
orem L. In fact, all will be special cases of the solutions constructed for arbitrary B 
by the Hermite-Krichever Ansatz. 

When £ > 2, many of the spectral polynomials will have degenerate roots if 
92, 93 € C are appropriately chosen. This means that, for example, a pair of the 2^+1 
band edge energies can be made to coincide by moving the modular parameter m G 
C\{0, 1} to one of a finite set of complex values. We indicate how to calculate these, 
or the corresponding values of Klein's absolute invariant J = 92/(92 ~ 27g|) G C. 

J is the more fundamental parameter, in algebraic geometry at least, since two 
elliptic curves are isomorphic (birationally equivalent) if and only if they have the 
same value of J. The m ^ J correspondence (12.3)1 maps C \ {0, 1} onto C, and 
it also maps m e M \ {0, 1} (in fact, m g (0, |]) onto J G [l,oo). Formally it is 
6-to-l. Each value of J corresponds to six values of m, with the exception of J = 
(i.e., 52 = 0), which corresponds to to = i ± ^i, and J ~ 1 (i.e., 93 — 0), which 
corresponds to m = — 1, i, 2. Elliptic curv es with J = 0, 1 a re called equianharmonic 
and lemniscatic, respectively ( Abramowit z fc Stegunll965[ §§ 18.13 and 18.15). Any 
equianharmonic curve has a triangular period lattice, with lu' /lu = c^^'^^/^ ^ and any 
lemniscatic curve has a square period lattice, with lu' /lo = ±i. 



(a) Lame polynomials 
The Lame polynomials are classified into species 1, 2, 3, 4 ( Whittaker fc WatsonI 



I1927L § 23.2). This is appropriate for some forms of the Lame equation, but for the 
elliptic-curve algebraic form, a more structured classification scheme is better. 

Definition. A solution of the Lame equation ()2.9)l on the elliptic curve -£92,93 is 
said to be a Lame polynomial of Type I if it is single-valued and of the form C{x) 
or D{x)y, where C, D are polynomials. A solution is said to be a Lame polynomial 
of Type II, associated with the branch point of the curve (7 = 1,2,3), if it is 
double-valued and of the form E{x)y/x — e-y or F{x)y / ^x — e^, where E,F are 
polynomials. The subtypes of Types I and II are species 1, 4 and 2, 3, respectively. 

To determine necessary conditions on i and B for there to be a nonzero Lame 
polynomial of each subtype, one may substitute the corresponding expression (C{x), 
etc.) into the Lame equation ()2.9)) . and work out a recurrence for the polynomial 
coefficients. This is similar to the approach of expanding in integer or half-integer 
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Table 2. Lame polynomials of Types I, II 

£ Type I solution: C{x; B , g2, ga) or D{x; B, g2, gs) y 

1 — 

2 x-lB 

6 

3 y 

4 x'-^Bx+{^^B'^i;g2) 
5 

6 x« _ ^Ba.^ + „ ^g,)^ + (-33k -B' + lilSffa - iffs) 

7 [a:'^±:Bx+j^,B'~l,g2]y 

8 ^4 _ _L.5^3 ^ (^_^52 _ ^32)^2 + (__^i33^+ -S_Bg2 - ^,gz)x 
+( 7413120 -B"^ ~ sTIso^^ga + iko^9i + ellga) 



Type 11 solution: E{x\ B,e^, g2, g3)^x - e~7 or F{x; B,e^,g2,gs) y/^/x- ^ 



a; 



2 y/^/^ 

3 ' 
4 
5 
6 
7 



' + (-52^ - \e^)x + (^B^ ^ ^^g^ _ ig2 _ y/y^ 
^.B + le^)a;2 + (-ij^B^ _ 1 fie^ + fe^ _ f g^)x 



\x- 



+ (~ 30-B ~ 2*^7")^ + (t560-B + eO-B'^T ~ 8*^7 ~ 10452)3; 

+ (-- lo^br-S' - aikg'e7 + + 3iio^g2 + ^e^ff2 - l^ffs)] y/ya; - 

powers of a; — e-y jWhittaker fc WatsorJll92'R § 23.41), though it leads to four-term 
rather than three-term recurrences. For the Type I solutions at least, the present ap- 
proach seems more natural, since they are not associated with any singular point e-y. 

If C{x) = E.cja;^ D{x) = Y.jdjX^, E{x) = e,a:^ and F{x) = Y.jf,^', 
substituting the expression for each species of solution into (|2.9(l and equating the 
coefficients of powers of x leads to the recurrence relations 

{2j-t){2j+l+l)cj-Bc,+i 

- U + 2)(j + 1)52 - (j + 2)(j + 3)53 Cj+3 = 0, 
{2j -£ + 3){2j + £ + 4) - B d^+i 

- [j + 2)(j + |).g2 - {j + 2)(j + 3).g3 d,+3 = 0, 
(2j - ^ + l)(2j + ^ + 2) + [(4j + 5)ey - B] e,+i 

+ [-{j + 1)52 + 4e2](j + 2) - + 2)(j + 3)53 6^+3 = 0, 
{2j -l + 2){2j+t + 3) /, + [-(4j + 7)ey - B] 

+ [-(j + 1)52 - 4e2](j + 2) /,+2 - {3 + 2)(j + 3)53 /,+3 = 0. 

It is easy to determine the integers £ for which C, Z?, E, F may be a polynomial 



(3.1) 
(3.2) 
(3.3) 
(3.4) 



Proposition 3.1. If £ > I is odd, nonzero Type I Lame polynomials of the fourth 
species and Type II ones of the second species can in principle be constructed from 
these recurrence relations, withdegD = {i—3)/2 anddegE = {£—l)/2 respectively. 
{The former assumes £ > 3.) If £ > 2 is even, nonzero Type I Lame polynomial of 
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Table 3. Lame spectral polynomials of Types I, II 
(Most of the ones with £ > 5 disagree with those pubUshed bv lGuerritor^ \l9Q^ .) 



£ L),(B;g2,gsj 



1 


1 




2 




- 3^2 


3 


B 




4 


B' 


- 52g2B + 560^3 


5 


B^ 


-27(72 


6 


B^ 


- 294^25^ + 777633-B + 3465^2 


7 


B^ 


- 196g2B + 2288^3 


8 


B^ 


- 1044^25^ + 48816535^ + 112320ff|B - 466560052^3 



£ Li {B;ej,g2,g3) 

1 B-e^ 

2 B + 3e^ 

3 B^- 66^,5 + (45e^-15p2) 

4 B^ + lOe-yB + (-35e^ - 752) 

5 B'^- 15e-,B^ + (315e2 - 13252)^ + (^e^52 + ^93) 

6 B'' + 21e-,B^ + (-189e2 - 8452)^ + (-^£^(72 + ^53) 

7 B* ~ 28e^B^ + (1134e^ - 57452)5^ + (3409e^p2 + 852533)5 

I / 292383 „2 „ 175175 „ „ _L00nQ„2\ 

8 + 366^5^ + (-594e^ - 414p2)B^ + (-9855e^p2 + 12285^3)5 

I ^ 245025 „2„ I 552825 „ _ i 7^oc;^2\ 

+1 4 e^g2+ 4 e^g3 + 7425ff2j 

the first species and Type II ones of the third species can be constructed similarly, 
with degC — £/2 and degF ^ {£— 2)/2 respectively. 

The coeiEcients in each Lame polynomial are computed from the appropriate 
recurrence relation by setting the coefficient of the highest power of x to unity, and 
working downward. Unless B is specially chosen, the coefficients of negative powers 
of X may be nonzero. But by examination, they will be zero if the coefficient of x~^ 
equals zero. 

Definition. The Type-I Lame spectral polynomial L\{B; §2, gs) is the polynomial 
monic in B which is proportional to the coefficient d-i ii£ is odd and c_i if ^ is even. 
(The former assumes ^ > 3; by convention L\ := 1.) The Type-II Lame spectral 
polynomial LY{B; e^, g2, gs) is similarly obtained from the coefficient e_i if £ is odd 
and /_i if ^ is even. Each spectral polynomial may be regarded as Yl^[B—Bs{g2,g3)], 
resp. Hsl-^ ~ Bs{e^, g2, gs)], where the roots {Bg} are the values of B for which a 
Lame equation solution of the indicated type exists, counted with multiplicity. 

By examination, iV^ :— degi^ is {£ — l)/2 if £ is odd and £/2 + 1 if ^ is even; 
and Nj^ := degL^^ is {£ + l)/2 if £ is odd and £/2 if £ is even. So as expected, 

3 

U{B;g2,g3) L\{B;g2,g3) J| L"(B; e^,, 52, 53), (3.5) 

7=1 

the full Lame spectral polynomial, has degree iV^ + 3iV^^ — 2£ + 1 in B. 

It should be noted that 117=1 -^^(-^5 ^7' 52, 33), the full Type-II Lame spec- 
tral polynomial, is a function only of B;g2,g3, since any symmetric polynomial 
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Table 4. Cohn polynomials of Types I, II 

£ Type I Cohn polynomial 

1 ~ 

2 J 

3 — 

4 2^ S'^ J + 5^ 7^ 

5 J 

6 2''3^5V^ + 11-37- 597-2^3^ 

7 2*3^5^ + 11^13^ 

8 J (2^2 3^^ 5^ 7^ + 2** 3^ 3664447 - 2'' 3^ - 397 - 364069 J + 113^^) 



£ Type II Cohn polynomial 

1 ^ 

2 — 

3 2^ J + 1 

4 2^ 3^^ J - 5^ 

5 2" 3"^ + 2^ 3^ 5^ 109 - 2^ 5* 17 - 151J + 5*^ 7^ 

6 2^^ 3*^ 5^ + 2" 3^ 17 - 359 + 2-* 57774169 J + 3^ 109^ 

7 2^° 3^^ 5'' + 2^"^ 3" S"* 19 - 22307 + 2" 3^^ 5^ 22158751 J* 

+ 2^ 3*^ 5^ 1276543 - 414016613 - 2^ 3^ 5^ 47202908378639011 
+ 3^" 29 - 41 - 101 - 895253 - 8050981 J - 2^ 5*^ 11^ 37^ 113^ 

8 2^° 3^ 5^ 7^ J® + 2^*^ 3^ 107 • 419 + 2^^ 3^ 12486499 J* - 2^ 3^" 1171 - 10477 

- 2-* 3'' 11 - 47 - 91938173 + 3* 20593 - 844499 J - 2^ 7^ 13^ 29^ 



in 61,62,63 can be written in terms of (72,33- For example, 616263 = 53/4. This is 
why 67 is absent on the left-hand side of H3.5|l . When using the recurrences (jH.lll - 
(I3.4|l . one should also note that 4ei^ — 9267 — 53 = 0, so = 1(5267+53). Any 
polynomial in 67,52,53 can be reduced to one which is of degree at most 2 in e^, 
much as any polynomial in x, y can be reduced to one of degree at most 1 in y. 

The Lame polynomials of Types I and II are listed in table El and the cor- 
responding spectral polynomials in table |21 They replace the table of fGuerrito3 
l|l909t ). with its many unfortunate errors. T he spectral polynoni ials with £ < 7 were 
recently computed by a different technique ijvan der Waall200l table A.3). The ta- 
ble of van der Waall displays the full Type II spectral polynomials, rather than the 
more fundamental 67-dependent polynomials LFf-B; e^, q^, q.^). The spectral polv - 
nomials can also be computed by the technique of iGesztesv & WeikardI l|l995oj) . 
which employs the Weierstrassian counterpa rt of the Ansatz used by He rmite in his 
solution of the Jacobi-form Lame equation ijWhittaker &: Watsonlll92'R §23.71). 

The roots of the spectral polynomials are the energies B for which the Lame 
polynomials are solutions of the Lame equation. It is clear that when i > 9, the 
Type II energies cannot be expressed in terms of radicals, since the degree of the 
spectral polynomial will be 5 or above. When £ = 8 or ^ > 10, the Type I energies 
cannot be so expressed. These statements apply also to the coefficients of the Lame 
polynomials, which depend on B. So when £ > 10, the symbolic computation of the 
Lame polynomials is impossible, and when £ = 8 or 9, it is possible only in part. But 
when 52 , 53 take on special values, what would otherwise be impossible may become 
possible. For instance, when 53 = (the lemniscatic case, including m = ^), the 
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quintic spectral polynomial L|(_B; 52,33) reduces to — 1044^25^ + 112320(7|i3, 
the roots of which can obviously be expressed in terms of radicals. 

In the context of the Jacobi form, the 2^+1 values {£^s(™)}slo ^'^^ which a Lame 
polynomial solution exists can be thought of as the 2£ + 1 branches of a spectral 
curve that lies over over the triply-p unctured s phere \ {0,1, 00}, the space of 
values of the modular parameter m. iTurbineJ l)l989l) showed that if £ > 2, this 
curve has only four connected components, not 2^+1. It is now clear why. These 
are the Type I component and the three Type II components, one associated with 
each point e-y. Since each of the four is defined by a polynomial in E and m, each 
can be extended to an algebraic curve ov er P-^ . At the v alues m = 0, 1, 00, the four 
curves may touch one another. (See, e.g.. lLi et alJ l)2000l fig ure 3), for the behavior 
of the real portions of the ^=1,2 curves as m 0, 1, and lAlhassid et all lll983l) 
for £ = 3.) These three values of m correspond to two of ei, 62, 63 coinciding, and 
the elliptic curve = 4:Ylj^i{x — e^) becoming rational rather than elliptic. Level 
crossings of this sort are perhaps less interesting than 'intra-curve' ones. 

In the present context, E is replaced by the transformed energy B, and m by the 
pair (72, 53 or the Klein invariant J, with J — 00 corresponding to m = 0, 1, 00. It is 
easy to determine which finite values of J yield coincident values of B. One simply 
computes the discriminants of the Type I and full Type II spectral polynomials, 
L\{B; g2, gs) and Oy^i ^7' 52, 53)- Each discriminant is zero if and only if 

there is a double root. By using J = 52/(52 ~ ^Tgl), one can eliminate 52,53 and 
obtain a polynomial equation for J. For each of Types I and II, there are coincident 
values of B if and only if J is a root of what may be called a Cohn polynomial. 

In table 01 the Cohn polynomials are listed. Since the coefficients are rather 
large integers that may have number-theoretic significance, each is given in a fully 
factored form. An interesting feature of these polynomials is that none has a zero 
on the real half-line [l,oo). Since J € [l,oo) corresponds to m € (0,1), the ex- 
istence of such a zero would imply that for some m g (0,1), two of the 2£ + 1 
band edge s become degenerate. That this cannot occur follows from a Sturmian 
argument llWhitt aker &: Watsonlll927L § 23.41). It also follows from the analysis of 
iGesztesv fc WeikardI lll995ol S 3V 

Proposition 3.2. For any integer £ > I, the degeneracies of the algebraic spec- 
trum of the Lame operator, which comprises the 2£ + I roots (up to multiplicity) of 
the spectral polynomial L^(_B; 52, 33), ^'"e fully captured by the Cohn polynomials of 
Types I and IL As the parameters 52,53 o-f£ varied, a pair of roots will coincide, 
reducing the number of distinct roots from 2£ + 1 to 2£, if and only if the Klein 
invariant J is a root of one of the two Cohn polynomials; and there are no multiple 
coincidences. 

This proposition will be proved in §0] The following conjecture is based on a close 
examination of the spectral and Cohn polynomials, for all £ < 25. 

Conjecture 3.3. 

L As a polynomial in J with integer coefficients, no Cohn polynomial has a 
nontrivial factor, except for the Type I Cohn polynomials with £ = 2 (mod 3), 
each of which is divisible by J. {These factors of J are visible in tabled) 
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2. If Ng and denote the degrees of the spectral polynomials of Types I, II, 
which are given above, then the Cohn polynomials of Types I, II have degrees 
{Nf - N} + 4)/6 I and NfiNf - l)/2, respectively. 



The conjectured degree formulas constitute a conjecture as to the number of points 
in elliptic moduli space (elliptic curve parameter space), labelled by J, at which the 
2£ + 1 distinct energies in the algebraic spectrum are reduced to 2^. For example, 
N\, equal 1, 2, so when £ ~ i the Cohn polynomials of Types 1,11 have degrees 
0, 1 respectively. That is, if £ = 3 there is no Type-I polynomial, and the Type-II one 
is linear in the invariant J. According to table 01 it equals 4J+ 1. A nondegeneracy 
co ndit ion equivalent to the linear condition 4 J + 1 was previously worked out 
bv lrTeibich 0994 §6.6), namely n^=i(5ff2 - 126^) ^ 0. 

The remarks regarding extra J factors amount to a conjecture that in the equian- 
harmonic case m — |± -^i (i.e., J = or 172 =0), there are only 2^ distinct energies 
if and only if £ = 2 (mod 3) . For those values of £, the double energy eigenvalue is 
evidently located at B ~ 0. It should be mentioned that when J = and £ = 
(mod 3), there is also an eigenvalue at B = 0, but it is a simple one. 

A periodicity of length 3 in ^ is present in the equianharmonic case of a third- 
order equation resembling H2.12|l , now called the Halphen equation llHalDhenlllSsi 
pp 571-4). By the preceding, a similar periodicity appears to be present in the 
equianharmonic case of the Lame equation. This was not previously realized. 







j> X exp 






. J y . 



(6) Twisted Lame polynomials 

The twisted Lame polynomials are exponentially modified Lame polynomials. 
They will play a major role in Theorem L and in the hyperelliptic reductions fol- 
lowing from the Hermite-Krichever Ansatz, but they are of independent interest. 

Definition. A solution of the Lame equation H2.9|l on the elliptic curve -Egj.sa is 
said to be a twisted Lame polynomial, of Type I or Type II associated with the 
point e^, if it has one of the two forms 

Cix) + Dix)y 
E{x)^x - + F{x) y/ 

with K G C nonzero. Here C, D, E, F are polynomials. 

On the level of Fuchsian differential equations, there is little to distinguish be- 
tween twisted Lame polynomials and ordinary Lame polynomials, which are simply 
twisted polynomials with k — 0. A function 'i'{x,y) — 5'(a;,?/)exp [k Jdx/y] will 
be a solution of the Lame equation (|2.9|l if and only if ^I* satisfies 

k2]|$ = o. (3.6) 

This is a Fuchsian equation on Eg^^g^ that generalizes but strongly resembles H2.9() . 
It has a single regular singular point, {x,y) = (00,00), and its characteristic ex- 
ponents there are {—£,£+ 1}. Like B, k is an accessory parameter that does not 
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affect the exponents. The values (B, k) for which a twisted or conventional Lame 
polynomial solution of H2.9() exists can be viewed as the points in a two-dimensional 
accessory parameter space at which has single or double- valued solutions. 

If C(x) = Dix) = Ejdjx^^ E{x) = and F{x) = fi^' ^ 

substituting the expression for each type of twisted polynomial solution into (|2.9|) 
and equating the coefficients of powers of x yields the coupled pairs of recurrences 

(2j - £)i2j +i+l)cj + {K^-B) c,+i 

- (j + 2)(i + 1)92 Cj+2 - (j + 2){j + 3)53 Cj+3 
+ 2k (4j -f- 2) - 2k {j + 1)52 rfj+i 

- 2k (j + 2)53 d,+2 = 0, 

(2j -£ + 3){2j + £ + 4) dj + {k^ - B) d^+i 

- (j + 2)(i + 1)52 dj+2 - {.] + 2)(j + 3)53 d,+3 
+ 2k {j + 2) c,+2 = 0; 

(2j - £ + 1) (2j + ^ + 2) + [(4i + 5)e^ + ^2 _ 5] e^.^^ 

+ [-{j + 1)92 + 4e2](j + 2) e,+2 - {] + 2)(j + 3)53 ej+3 
+ 2k (4j + 4) /, -f 2k (4j + 6)e^ 
+ 2.^ (4j + 8)(e2-i52)/,+2 = 0, 

(2i -l + 2){2j+£ + 3) + [-(4j + 7)e^ + - i?] 

+ [-(j + i).92 - 4e2](j + 2) /,+2 - (j + 2)(j + 3)53 /,+3 
-I- 2k (j 4- I) gj+i - 2k (j 2)e-y ej+2 = 0. 

It is easy to determine the maximum value of the exponent j in each of C, D, E, F. 

Proposition 3.4. Nonzero twisted Lame polynomials of Type I {if £ > 3) and of 

Type II {if £ > 2) can in principle be constructed from these recurrence relations. If 
£ is odd, resp. even, then degC\ D; E, F are {£~l)/2,{£-3)/2;{£-l)/2,{£-3)/2, 
resp. £/2, {£ — A)/2\£/2 — l,£/2 — 1. The coefficients are computed by setting the 
coefficient of the highest power of x to unity in D and E, resp. C and F, and 
working downward. 

Unless B, k are specially chosen, the coefhcients of negative powers of x may 
be nonzero. But by examination, they will be zero if the coefficients of x^^ in C 
and D (for Type I), or E and F (for Type II), equal zero. c_i = 0, d_i — 0, and 
e_i = 0, /_i ~ 0, are coupled polynomial equations in B,k, and their solutions 
may be computed by polynomial elimination, e.g., by computing resultants. A minor 
problem is the proper handling of the case k = 0, in which (|3. 7(1 ~ (|3. 10(1 reduce to 
((3. 1(1 - 1(3.4(1 . If £ is odd, resp. even, then c-i and /-i, resp. d-i and e-i, turn out to 
be divisible by k. By dividing the appropriate equations by k before solving each 
pair of coupled equations, the spurious k = solutions can be eliminated. 

Definition. The Type-I twisted Lame spectral polynomial Lt\{B; 92, 93) is the 
polynomial monic in B which is proportional to the resultant of c-i,d-i with re- 
spect to K, with K factors removed as indicated. (This assumes ^ > 3; by convention 
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Lt\ = Lt\ :— 1.) The Type-II twisted Lame spectral polynomial LtY{B; e^, 172,53) 
is similarly obtained from e_i, /-i. (This assumes € > 2; by convention Lt^ ■= 1-) 
Each twisted spectral polynomial may be regarded as Y\g[B — i3s(ff2,33)], resp. 

— Bs{e^, g2, 53)], where the roots {Bg} are the values of B for which a Lame 
equation solution of the specified type exists, counted with multiplicity. 

The twisted Lame spectral polynomials for £ < 8 are listed in table El The 
polynomials Lt\^ and Lt^^ are omitted on account of lack of space (their respective 
degrees are 12 and 16). The following proposition will be proved in §0] 

Proposition 3.5. 

1. For any integer £ > 3, resp. i > 2, there is a nontrivial twisted spectral 
polynomial 0/ Type I, resp. Type IL 

2. For any integer £>1, Nt\ := deg Lt\ is {£'^ - l)/4 if £ is odd and i'^/A - 1 if 
£ is even; and NtY '■= degL^^ is (€^ — l)/4 if £ is odd and £^/4 if £ is even. 

So the full twisted Lame spectral polynomial Lte{B] g2, gs), which by definition 
equals Lt\{B; 52, 33) D^^i Ltf{B; e-,, 52, 53), wiU be of degree Nt\+2,Ntf = £"^-1 in 
the spectral parameter B. Like the ordinary degree-(2£+ 1) spectral polynomial L^, 
Ltf is a function of _B; 52,53 only, because any symmetric polynomial in 61,62,63 
can be written in terms of the invariants (72,53- 

(c) Theta-twisted Lame polynomials 

Lame equation solutions of a third sort can be constructed for certain val- 
ues of the spectral parameter B. These are linear combinations, over polynomials 
in the coordinate x, of (i) the multi-valued meromorphic function $(a;, y; Xq, j/o) 
parametrized by the point (a;o,2/o) G -^92.33 \ {(00,0°)}, ^-i^d (ii) its derivative 

$(^)(a;,y;xo,2/o) := {v^^ ^{x,y;xo,ya) = ^ {j^rZ~^ ^{x,y; XQ^yo). (3.11) 

One way of seeing that $, <i>(^) are a natural basis is to note that when {xq, yo) = 
{ej,0), they reduce to ^/x — e-y, \yl^Jx — e^. So the class of functions constructed 
from them will include the Lame polynomials of Type 11. 

Definition. A solution of the Lame equation (|2.9|l on the elliptic curve -Ega.ss is 
said to be a theta-twisted Lame polynomial if it is of the form A{x)(^{x, y; xq, yo) + 
2B{x)^^^'>{x,y;xo,yo), with {xo,yo) 7^ (6^,0) for 7 = 1,2,3. Here A,B are polyno- 
mials, and the innocuous '2' factor compensates for the 'i' factor of (|3.11|) . 

If A{x) = ttjX^ and B{x) = hjX^ , substituting this expression into (|2.9|l 
and equating the coefficients of powers of x yields the coupled pair of recurrences 

(2j ~£+l){2j + £ + 2) a, + [(4j + 5)a:o - B] a^+i 

+ [-{] + 1)92 + 4x2] (j ^ 2) a,+2 - {] + 2)ij + 3)53 a,+3 (3.12) 
- 22/0 (4j + 6) bj+i - 4x02/0 (j + 2) = 0, 

(2j -£ + 2){2j+£ + 3) b, + [-{Aj + 7).xo - B] b,+i 

+ [-(.? + 1)52 - 4xg](j + 2) 6,+2 - (.7 + 2)ij + 3)53 b,+3 (3.13) 
+ yoU + 2) aj+2 = 0. 
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Table 5. Twisted Lame spectral polynomials of Types I, II 



Lt\{B;g2,g3) 



1 
1 

+ ( 

+ 
+ 

+ 
+ 
+ 
+ 
+ 
+ 

+ 
+ 
+ 



897 „ o4 19845 „ d3 , 546993 „2 p! 
9?; ."^ ?;nR:?4?;i 7f; " 



2 2 
4100625 „3 506345175 „2n 

, 4 52 16 93) 

2751 „ d6 181521 „ r>5 j 

—2-92B 2 93B -\ Jg 

' " 16273476669 2-102 
16 ' 



2 T-,2 _|_ 



3407481 „2 p4 , 
16 g2-^ + 



™e5si-— 16--53;- -- 

16980286476 „4 , 1664232587326 2\ 

4 92 + 16 5253) 



4186ff2Sl° -'i04|223g3B9 ^ 171^3^2^. ^ 

4590448625 „3 59437238487 „2\ p6 3146848477773 „2 „ pi 
4 92 2 5'3)'° 32 ff253-0 



/2-" -1 8 h 

3362086035 g2g^£ 

3510785355 g^g3g7 



4 52 2 9S)^ 32 

239496271862939 „4 , 349377693363699 „ „2-, p4 

256 92 + 16 52S3)S 

167403005460693 „3 „ 10531741687878125 „3^ p3 
77^ 9-2 93 — 53 )^ 



256 92 i 16 i 

2167403005460693 „3 „ 105317416878781- , „ -- 

128 S'2S'3 32 9a )^ 

1196552376313749 „S 243386984019562383 „2,.2', p2 

8 92 64 92g?,)B 

570084251356448829 „4„ , 15458554942852896875 „ „3\ n 
128 92 93 H 128 92g3)B 

1649721227262688125 „6 , 16766233150463677881 „3„2 , 285799721595172159375 „4\ 

256 92 H 128 5253 H 256 93) 

OI8852S" _ 4944861 ^^512 ^ 48623733 ^2^11 ^ 33098210361 ^^^^^10 
r.r.,,,^^,.„ , ,,-..„ ,)8i9345i 2-, p9 46667883177495 „2 „ p8 
53)^ 32 5253-D 



^^^^^y^^ — — ysu -I 

210211163145 „3 1634908193451 ,.2\ p,. . „„„ 

, 8 92 4 53ji^ 32 

• 9879747455405475 „4 , 13114846350610875 „ „2\ p' 

. 266 92 + 16 5253)'« 

■ 2796144976900437S „3 „ 268903662388069375 „3^ p6 
, 128 5253 32 53j« 



■g^9i)i 



>B^ 



^ 266 92 H 16 52 

r 2796144976900437S 3 _ 268903662388069375 

k 128 52 53 32 

r 1268095592996251875 „5 160577779706139(i5(i2 

^ 92 

r 6731()i«21U«52!)i«4:i75 ,.4„ , 19<)Si2(i47."')8rjr)(i!)!)190625 „ „3\ p4 
. 128 f253 H 128 5253)« 

f 1247302375822866515625 „6 , 15203913100824300328125 „3 „2 , 83437068242769811171876 „4^ p3 

256 92 H 128 5253 H 256 93)^ 

75 -2 -3\ p2 



128 

124y3U23yb822S6651b625 „0 , 15203913100824300328125 „3 „2 I »ci43YUt>»:J4iiYt>y»lilYl»Y0 „4\ pd 

256 92 H 128 5253 H 256 93)^ 

2476060022819411015625 „5 „ 67800902314274734921875 „2„3\p2 
Ifi 5253 16 9293)B 

8803296317899887890625 „7 100408533824S6687S390625 „4 „2 1003123416299777251171875 „ „4\ p 
16 52 8 5253 16 92g3)B 

14047813273244501953125 gig^ _ 305679140836374550781250t?|ffi 



— Ifi 52y3 o.^o.- 

8803296317899887890625 „7 100408533824S6587S390625 „4 „2 101 
16 52 8 5253 

,7813273244501953125 gig^ _ 305679140836374550781250t?; 

6263785284763018974609375 „5^ 

2 —53; 



Lt'e{B;e^,g2,g3) 



B — Qe~f 

B'^ - 15e^B + (-225e2 

B* - 55e^B^ + (-945e^ + ^52)^^ + (1960e^52 + 2450^3)5 
+ (61740e^32 



75 \ 
—52) 
539 , 
4 



■ 686006^53 - 9261g'i) 



+ ( 



lOSe^B 

6077485 



' + (- 
6^52 - 



■7245e^ + 



lYUY „ \ r4 , /] 

-^92)B +{- 



-e-y52 + 



'-93)B' 



-e-ygs 



16 



^52^)5^ 



4 ^7 _ 
I / 56260575 ^2 „ _i_ 1 CQnonrio „2 I 64117316 „ „ 

+ ( 4 e^53 + 153U9UUe7ff2 H tr — 5253 



4 T 
I / 86113125 „2„2 
+ I 2 '^752 



- 231e. 



2 



+ (- 

,2 , 



72454095 „ 

4 ^752 - 

118721295 „2 „ 
6^53 



2 

21735e^ + 

58128273g 



16 

5„ „ „ I 36905625 „3 , 

-675253 H 8 52 + 

6699 „ \ r7 I / 255087 

-^92>B 



506345175 „2 



gl) 



4 

7116417 



53 

gl 



+ (— 

38437119^2^1^5 



6752 + 



-53 



]B'^ 



16 

249847983 



5253 



+ ( 4 ^-,y2 4 '^7. 

+ ( 4 '-7; 

, f 8838982845 ^2 „2 . 10241858781 

+ \ 2 ^^52 H 2 ^7 

+ ( 6418.527885 g2g^^^ ^ 8690^5655 ^^^3 

_^ ( 1189804891275 ^,2 ^3 ^ 1089315875340e^ 1906302781845 



675253 + 



4499132715 ^3 ^ 2338773426fli ) 



-8637504668167 33 - 



;53 



13550225535 „2 



783(13398827.1 



e-ygigs 



5253)S" 



16 

, / 11649628111276 „2 „2„ 
+ ( 4 e^5253 



c/2 - 1512938715755253)5 



33558599497567 ffl + 499269776197567^253 



16101369773875 ^3^^ _^ 128383656736503| ) 
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Table 6. Theta-twisted Lame spectral polynomials 



L0i{B-g2,g3) 



1 1 

2 1 

3 1 

5 - Iflg^B - 2^ff3 

6 B"^ - ^g2B^ - i2fS5g3B3 + leooSSffiS^ + ^^^^^gigsB 

+ (- 57635^606 ^2 ^ 96850215ff|) 

7 58 _ 19|65^^56 _ 83|43g3B5 + 26047931 ^2^4 ^ 4205970769 ^^^^5 3 

I 1 37048456991 „3 204966441251 „2'i p2 , 8684628953 „2 „ o 

+ ( 48 f2 16 53J-D H 6 9293B 

, I 552623218875 „4 , 43902444356771 „ J2\ 

+ ( 4 92 H 12 92^3) 

8 5I2 _ 18063^2^10 _ 4067739^^^9 ^ 73174185^2^8 ^ 22697632971 ^^^^^7 

+ 2^ ( 102471 ISffl - 167402573<7i)B« - 1385229823965 ^2^^^5 

+ 16422883331 + 360649430752^3)5* 
+ 63969750000 (2175^153 - 101062g|)B^ 
+ 98415000 (62738863gf - 165603184531 
+ 921164400000 (-256O363I33 + 7098507g2ff3)B 
+ 19683000000 (-35153041g| + 669725199ffigg + 7578832716gt) 



Proposition 3.6. If i > nonzero theta-twisted polynomials can in principle be 
computed from these recurrences. If I is odd, resp. even, then degA,B are [l — l)/2, 
{£—5)/2, resp. £/2—2,i/2^1. The coefficients are computed by setting the coefficient 
of the highest power of x in A, resp. B, to unity, and working downward. 

Unless B and the point (xq, j/o) ^^re specially chosen, the coefficients of negative 
powers of x may be nonzero. But by examination, they will be zero if the coeffi- 
cients of x~^ in A and B are both zero. a_i = 0, 5_i = are equations in B; xq, y^. 
Together with the identity j/q = 4a;o — 1722^0 — 53 j they make up a set of three polyno- 
mial equations for these three unknowns. This system may be solved by polynomial 
elimination. For example, to obtain a single polynomial equation for B (involv- 
ing g2,ff3 of course), one may eliminate j/o from a_i = 0, 5_i = by computing 
their resultants against the third equation; a nd then eliminat e xq. Alternatively, a 
Grobner basis calculation may be performed l|Brezhnevll2"ool . 

Irrespective of which procedure is followed, there is a minor problem: the han- 
dling of the improper case (xq, yo) — (e-y, 0), in which H3.12|l ~ H3.i;^(l reduce to H3.8|l - 
(|3.4|1 . If t is odd, resp. even, then the left-hand side of the equation 6_i =0, resp. 
a_i — 0, turns out to be divisible by y^. By dividing the appropriate equation by j/o 
before eliminating xo,yo, the spurious solutions with yo ~ can be eliminated. 

Definition. The theta-twisted Lame spectral polynomial LOi{B\ 52, ffa) is the poly- 
nomial monic in B which is obtained by eliminating xo,yo from the equations 
a_i = 0, 6_i = 0, with yo factors removed as indicated. (This assumes ^ > 4; by 
convention L9i — L62 — L9^ :— 1.) Each theta-twisted spectral polynomial may be 
regarded as Osl^ ~ ^s(52, ffa)], where the roots {Bg] are the values of B for which 
a theta-twisted Lame polynomial exists, counted with multiplicity. 

The theta-twisted Lame spectral polynomials for ^ < 8 are listed in table El The 
following proposition will be proved in §0] 

Proposition 3.7. 
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L For any integer i > A there is a nontrivial theta-twisted spectral polyno- 
mial LOi. 



2. For any integer i > 2, NOg := degLOi is 3)/4 if £ is odd and 

i{£-2)/A ife is even. 



4. The Hermite— Krichever Ansatz 



The Hermite-Krichever Ansatz is a tool for solving any Schrodinger-like differential 
equation, not necessarily of second order, with coefficient functions that are elliptic. 
Such an equation should ideally have one or more independent solutions which, 
according to the Ansatz, are expressible as finite series in the derivatives of an 
elliptic Baker-Akhiezer function, including an exponential factor. (Cf. (ILSfl .) 

In the context of (|2.9|l . the elliptic-curve algebraic form of the Lame equation, 
this means that one hopes to be able to construct a solution on the curve Eg^ g^, 
except at a finite number of values of the spectral parameter i? G C, as a finite series 
in the functions ^^^'> :— (yd/da;)-' $, j > 0, multiplied by a factor exp [n J da;/y]. 
Here <&(•,•; xq, j/o) is the fundamental multi- valued meromorphic function on Eg^_g.^ 
introduced in §|5] Actually, a different but equivalent sort of series is easier to ma- 
nipulate symbolically. By examination <i>^^^ — {2x + a;o)$, from which it follows by 
induction on j that any finite series in j > 0, is a combination (over polyno- 
mials in x) of the basis functions <&, <i>(^^. This motivates the following definition. 



Definition. A solution of the Lame equation (|2.9|) on the elliptic curve Eg^^ g^ is 
said to be an Hermite-Krichever solution if it is of the form 



Aix)<S>{x, y; XQ,yo) + 2B{x)^^^Hx, y; xq, yo) 

y + ya 



A{x) + B{x) 



X — Xq 



exp 

^{x,y;xo,yo) exp 



dx 

y _ 

dx 

y _ 



(4.1) 



for some {xq, yo) e Eg^^^g^ \ {(oo, oo)} and k G C. Here A, B are polynomials. 

As defined, Hermite-Krichever solutions subsume most of the solutions explored 
in §13 If K = 0, they reduce to theta-twisted Lame polynomials. If (a^oi 2/o) = (67, 0) 
for 7 = 1,2,3, in which case $,<I>(^) degenerate to y/x — e-y, \y/ y/x — e-y, they 
reduce to twisted Lame polynomials of Type II. If both specializations are applied, 
they reduce to ordinary Lame polynomials of Type II. The Lame polynomials of 
Type I, both ordinary and twisted, are not of the Hermite-Krichever form, but they 
can be viewed as arising from a passage to the (xq, yo) — > (00, 00) limit. 

If A{x) = J2j O'j^'' and B{x) ~ J2j bjX^ , substituting (|4.1|) into 12. 9|) and equat- 
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ing the coefficients of powers of x yields the coupled pair of recurrences 

{2j -£+l){2j + e + 2)aj + [(4j + 5)xo + ~ B] a^+i 

+ [-{j + 1)92 + - 2Kyo]{] + 2) a,+2 - (j + 2)(j + 8)33 a^+s 
+ 8k{j + 1) bj + A{kxq - yo){2j + 3) 
+ 2 [/t(4a;2 - g^) - 2xoyo] {j + 2) &j+2 = 0, 
(2j -£ + 2)i2j+i + 3) fe, + [-(4j + 7)xo + - B] 

+ [-{j + 1)32 - 4a;2 + 2Kyo]{3 + 2) 6,+2 - {j + 2)(j + 3)53 6j+3 (4-3) 
+ K(2j + 3) flj+i - (2kxo - yo){j + 2) 0^+2 = 0. 

If K 0, reduce to (|TT^ - (|^T^ . and if (xcyo) = (£7,0), they reduce 

to (|3.9|l - (|3.10() . If both specializations are applied, they reduce to (|3.3|I - H3.4|I . 

Proposition 4.1. For all i > 2, Hermite-Krichever solutions can in principle be 
computed from these recurrences. If I is odd, resp. even, then deg A, B are {£ — l)/2, 
(£— 3)/2, resp. i?/2— 1, £/2— 1. The coefficients are computed by setting the coefficient 
of the highest power of x in A, resp. B, to unity, and working downward. 

Unless B, k and the point (xQ^yo) are specially chosen, the coefficients of neg- 
ative powers of x may be nonzero. But by examination, they will be zero if the 
coefficients of x~^ in A and B are both zero. a_i = 0, 6_i = are equations 
in B] K;xo,yo. They are 'compatibility conditions' similar to those that appear 
in other applications of the Hermite-Krichevcr Ansatz. Together with the identity 
?/q = 4:Xq — g2Xo — gs, they make up a set of three equations for these four unknowns. 

Informally, one may eliminate any two of B; n; xo,yo, and derive an algebraic 
relation between the remaining two unknowns (involving 52, 33 of course). A rigor- 
ous investigation must be more careful. For example, if the ideal generated by the 
three equations contained a polynomial involving only B (and 32,33), th en a solu- 
tion o f the Hermite-Krichever form would exist for very few values of B ((BrezhnevI 
In practice, this problem does not arise : except at a finite number of values 
of B (at most), the Ansatz can be employed ijGesztesv fc Weikardlll998| . In fact, 
in previous work an algebraic curve in [B, k) has been derived for each £ < 5. The 
{B, K)-curve is the one with the most physical significance, since i? is a transformed 
energy and n is related to the crystal momentum. 

Solving for (a;o,yo) as functions of {B,k) reveals that xq is a rational function 
of B, and that if k is not identically zero, j/o is a rational function of B, times k. 
These facts can be interpreted in terms of the following seemingly different curve. 

Definition. The ^'th Lame spectral curve := F^ ((72, 33) is the hyperelliptic curve 
over 3 B comprising all {B,v) satisfying = L£(_B; 52, 33), where is the full 
Lame spectral polynomial, of degree 2£ -\- 1 in B. {T( was informally introduced 
as F^ in §n where the original energy parameter E was used.) F^ will have genus £ 
unless two roots of L£(-; 52,33) coincide, i.e., unless the Klein invariant J is a root 
of one of the two Cohn polynomials of table 21 in which case the genus equals £—1. 

For each £, there must exist a parametrization of Lame equation solutions by a 
point [B, v) on the punctured curve Te \ {(oo.o o)}. by the general theory of Hill's 
equation on R ((McKean & van MoerbekgllQTfll) . For any finite-band Schrodinger 
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equation on an elliptic curve, including the integer-^ Lame equation, the Baker- 
Akhiezer function (|2.13|) provides such a parametrization of solutions. In the general 
theory, the parametrizing hyperelliptic curve T for any finite-band Hill's equation 
arises from differential-difference bispectrality: as a uniformization of the relation 
between the energy parameter B and the crystal momentum k (Trcibich 2001). This 
curve r 9 {B, v) is defined by an irrationality of the form = \-{B), and B and k 
are meromorphic functions on it; the former single- valued, and the latter additively 
multi-valued. The energy is computed from the degree-2 map _B : F — > given by 
{B^v) _B, and the crystal momentum from the formula 

--'/[Kf^)-<^-']f//f' 

in which the line integrals on Eg^^g^ are taken over the appropriate fundamental 
loop. Here {BjV) i-^ {xo,yo) is a certain projection n : T ^ ^92,93^ ^^'^ k : T ^ 
is a certain auxiliary meromorphic function. These two morphisms of complex man- 
ifolds are 'odd' under the involution {B,u) i-^ {B,—v), i.e., xoiB^—v) = xo{B,i') 
and yo{B,—u) = —yo{B,u), and k{B,—v) — —k{B,u). So xq must be a rational 
function of B, and each of yo and k must be a rational fimction of B, times v. 

In the general theory of finite-band equations, the Baker-Akhiezer uniformiza- 
tion is viewed as more fundamental than the Hermite-Krichever Ansatz. However, 
in the case of the integer-^ Lame equation one can immediately identify the curve 
in {B, k), derived from the Ansatz as explained above, with the ^'th Lame spectral 
curve 3 {B, v). It is isomorphic to it by a birational equivalence of a simple kind: 
the ratio k/v is a rational function of B. 

This interpretation makes possible a geometrical understanding of each of the 
types of Lame spectral polynomial worked out in §13 Due to oddness, each of the 
finite Weierstrass points {(Bs,0)}^£o on T^, which correspond to band edges, must 
be mapped by the projection tt^ : ^ ^92,93 to one of the finite Weierstrass points 
{(e-y, 0)}i^^j or to (00,00). Bearing in mind that Type I Lame polynomials, both 
ordinary and twisted, are not of the Hermite-Krichever form, on account of (xq, yo) 
formally equalling (00,00) and k equalling 00, one has the following proposition. 

Proposition 4.2. 

1. The roots of the Type- 1 Lame spectral polynomial L\{B) are the B-values of the 
finite Weierstrass points {{Bs,0)}'^^q that are projected by tt^ : F^ — + ^92, 93 
to the infinite Weierstrass point (00,00). Moreover, the roots of the Type-I 
twisted Lame spectral polynomial Lt\{B) include the B-values of the finite 
non- Weierstrass points that are projected to (00,00), and do not include the 
B-value of any finite point that is not projected to (00, 00). 

2. For 7 = 1,2,3, the roots of the Type-II Lame spectral polynomial LY{B]e^) 
are the B-values of the finite Weierstrass points { [Bg , 0)}slo ^^"^^ ^'^s projected 
by TTf : Ff — > Eg^^g^ to the finite Weierstrass point (e-y,0). Moreover, the roots 
of the Type-II twisted Lame spectral polynomial Lt^iB; Cj) include the B- 
values of the finite non- Weierstrass points that are projected to (e^,0), and 
do not include the B-value of any finite point that is not projected to (e-y, 0). 
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3. The roots of the theta-twisted Lame spectral polynomial L9i(B) include the 
B -values of the finite non-Weierstrass points that are zeroes of ni : Tg , 
and do not include the B-value of any finite point that is not a zero. 

Remark. The phrasing of the proposition leaves open the possibihties that (1) a 
root of Lt\{B) may be a root of L\{B), (2) a root of LtY{B,e^) may be a root 
of LY{B,e^), and (3) a root of L9\{B) may be the i3-value of a finite Weierstrass 
point, i.e., a band edge. Generically these three types of coincidence do not occur, 
but instances are not difficult to find. One is the case £ = Q (mod 3) and 32 = 
(i.e., J = 0), in which by examination Lt\ and L\ have the common root B = 0. 

To proceed beyond the proposition, a significant result from finite-band integra- 
tion theory is needed. For the integer-£ Lame equation, the covering ni : Ti -^32-33 
and the auxiliary map Kg :Ti ^ ¥^ are both of degree £{£ + l)/2, irrespective of the 
choice of elliptic curve Eg^^g^. From this fact, supplemented by proposition 14. 21 the 
two propositions 13.51 and 13.71 which were left unproved in §13 immediately follow. 

The fibre over any point {xq, j/o) G Eg^^g.^ must comprise £{£+ l)/2 points of Ti, 
the counting being up to multiplicity. Consider the fibre over (00,00), which by 
examination includes with unit multiplicity the point {B,v) = (00,00). It also 
includes each finite Weierstrass point of Vg that corresponds to a Type-I Lame 
polynomial. As was shown in §13 the number of these up to multiplicity, N\ :~ 
degL^, equals (£ — l)/2 if £ is odd and £/2 -I- 1 if £ is even. So the number of 
additional points above (00,00) \s £{£+!) /2- 1 -{£- I) /2 = -1)12 if£ is odd 
and £(£+\)l2-\~ {£/2 + l) = £'^/2-2 if £ is even. Since the projection tt^ is odd 
under the involution (i?, v) [B, — i^), these occur in pairs. So Nt\ := degLt^ must 
equal (£^ — l)/4 if £ is odd and £^/4 — 1 if £ is even; which is the formula for Nt\ 
stated in proposition 13.51 A similar computation applied to the fiber above any 
finite Weierstrass point (e-y, 0) yields the formula for Nt^ given in that proposition. 

The formula for N9g :— deg LOg stated in proposition 13.71 can also be derived 
with the aid of proposition^^! Since the map k : ^ has degree £{£-\-l)/2, the 
fibre above comprises that number of points, up to multiplicity. It includes each 
finite Weierstrass point of F^ that corresponds to a Type-II Lame polynomial (but 
not the Weierstrass points corresponding to Type-I Lame polynomials, since those 
are not of the Hermite-Krichever form and formally have k — 00). The number 
of these up to multiplicity is three times ■= degL^, which equals 3{£ + l)/2 
if £ is odd and 3^/2 if £ is even. So the number of additional points above is 
£{£+l)/2-3{£ + l)/2 = {£+l){£-3)/2 if £ is odd and £(£-M)/2-3£/2 = £{£-2)/2 
if £ is even. They conic in pairs, and division by two yields for the formula for NOg 
given in proposition 13. 71 

A geometrized version of proposition 13 . 21 can also be proved, with the aid of an 
additional result that goes beyond the Hermite-Krichever Ansatz. The Lame spec- 
tral curve F;((72,ff3) *s nonsingular with genus £ for generic values of the Klein in- 
variant J — ^(52,53)7 0''r>'d when it degenerates to a singular curve F;° := Ff((72,33), 
the singular curve has genus £ — I, with singularities that are limits as ((72,53) 
(52,53) 0/ Weierstrass points of Tg{g2,g3)- It follows that the £'th Type-I and 
Type-II Cohn polynomials, which characterize the pairs (52,53) for which the tth 
Lame operator on Eg^^g^ has degenerate algebraic spectrum of the specified type, 
i.e., for which F^((72 , 53) has a pair of degenerate finite Weierstrass points of the 
specified type, in fact do more: they characterize the (52,53) for which Tg{g2,g3) 
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is singular. Due to the reduction of the genus by at most unity, there can be, as 
proposition 13 . 21 states, no multiple coincidences of the algebraic spectrum. 

The problem of explicitly constructing an Hermite-Krichever solution of the 
integer-^ Lame equation, of the form (|4.1|) . will now be considered. What are needed 
are the quantities XQ,yo,K, or equivalcntly Xo,yQ/iy,K/v. Each of the latter is a 
rational function of the spectral parameter B. 

One way of deriving these functions is to eliminate variables from the system 
of three polynomial equations in B; k; xq, yo, as explained above. Coupled with the 
spectral equation = L^(_B), this yields exphcit expressions for the three desired 
functions. By the standards of polynomial elimination algorithms, this procedure 
is not time-consuming. It is much more efficient than the manipulation s of com- 
patibi lity conditions that other authors have employed. Beginning with iHalphenl 
l|1888|) . it has been the universal practice to apply the Hermite-Krichever Ansatz 
to the Weierstrassian form of the Lame equation, i.e., to (|2.()|l . rather than to the 
elliptic curve form 12.9|l . It is now clear that this Weierstrassian approach is far 
from optimal. For example, the computation of the covering map (B, v) i— > (xq, yo) 
in the case ^ = 5, by Eilbeck & Enol'skii (199'3), required seven hours of computer 
time. The just-sketched elliptic curve approach requires only a fraction of a second. 

It turns out that for constructing Hermite-Krichever solutions, this revised elim- 
ination scheme is also not optimal. Remarkably, at this point no elimination needs 
to he performed at all, since the covering tt^ : — > Eg^^g^ and auxiliary function 



Ke ■ Ff 



can be computed directly from the spectral polynomials of §|31 



Theorem L. For all integer £ > 1, the covering map tt^ : F^(g2,53) ^92,33 
appearing in the Ansatz maps {B,v) to (xa^yo) according to 



2^0(5; 52, 53) = -f 
2/o(B,i/;52,53) 



Lf{B; e^, g2, 33) [Ltf{B; e^, 52, ffs)] 



W+IW L\iB;g2,g3)[Lt\iB;g2,g3W 
16 / Ul=iLtf{B;e^,g2,g3) \ 



m+lW [L\{B;g2,g3W[Lt\{B;g2,g3W 



(4.5) 
(4.6) 



with 7 in H4.5|l being any of 1, 2, 3. The auxiliary function Ke : Te ^ is given by 



K{B,i^;g2,g3) 



L0i{B;g2,g3) 



L\{B;g2,g3)Lt\{B]g2,g3) 



(4.7) 



Proof. With the exception of the three i!-dependent prefactors, such as A/[£{£+l)Y , 
the formulas (|4.5|l - 14.7|l follow uniquely from proposition 14.21 regarded as a list of 
properties that tti and must satisfy. 

TT^ must map each point (-B^, 0), where Bg is a root of L\{B), singly to (00, 00), 
and each point (i3s,0), where Bs is a root of Lf{B] e-y), singly to (e^,0). It must 
also map each point (Btjiiyt), where Bt is a root of Lt\[B), singly to (00,00), and 
each point {Bt, ii^t), where Bt is a root of LtY{B; e^), singly to [e^, 0). In all these 
statements, the counting is up to multiplicity. 

Also, B ^ k/v must map each point {B',±v'), where B' is a root of L9g{B), 
singly to zero, and must map each point {Bg, 0), where Bg is a root of L\{B), and 
each point (Bt, ±i/(), where Bt is a root of Lt\{B), singly to 00. In these statements 
as well, the counting is up to multiplicity. 
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The i?-dependent prefactors in H4.5|) - (|4.7|) can be deduced from the leading-order 
asymptotic behavior of Xq and k/v as B oo. □ 

The remarkably simple formulas of the theorem permit Hcrmite-Krichever solu- 
tions of the form (|4.1|) to be constructed for quite large values of £, since the Lame 
spectral polynomials L, Lt, L6 (ordinary, twisted and theta-twisted) are relatively 
easy to work out, as §|3|made clear. Tables O O and |S1 may be consulted. 

It should be stressed that in the formula (|4.5|1 for xq, the same right-hand side 
results, irrespective of which of the three values of 7 is chosen. All terms explicitly 
involving e-y will cancel. Of course, all powers of higher than the second must 
first be rewritten in terms of 52,53 by using the identity = j{g2£j + 53)- In the 
same way, it is understood that the numerator of the right-hand side of (|4.6|l . the 
terms of which are symmetric in 61,62,63, should be rewritten in terms of 52,53- 
(This can always be done: for example, 6^62 + e\e\ -|- e|ef equals 5I/I6.) 

The application of Theorem L to the cases ^ = 1, 2, 3 may be illuminating. 

o If ^ = 1, then [xQ^yo) — {B,2v) and k = 0. The map tti : Fi ^ ^92.93 ^ 
mere change of normalization, since Fi is isomorphic to Eg^^g^] cf. H2.15|l . 

o lie = 2, then 

1 iB + 3e,)iB-6e, f B^ + 2753 

^''='^^+9 =9(52-352)' ^^-^^ 

_ 2 n^^i(^-6e^) _ 2(^3 _ Qg^B _ 54g^) 
^° 27 (52-352)2 " 27(52-352)2 ^ ' 

and K = - {2 / [3(52 - 352)]} z^. 

o If ^ = 3, then 



x„ 



l_ (52 - 66^5 + 45e2 - 1552)(5^ - 15e^5 - 2256.2 + ^g^) 
36 B{B^ - f 52)2 



165^ + 36O5254 + 27000535^ - 33755^5 

'"3 



3O375O52535 - 843755I -I- 22781255:2 



36 5(452 _ 75^2)2 

1 n'=i(^'-156^^" 225e2H-fg2) 
^° 108 B^{B^-^g2)^ " 

I656 - 1800525-* - 54OOO5353 

- I68755I52 + 42187551 - 113906255? 



(4.10) 



(4.11) 



2752(452-7552)3 
and K = - { 10 / [35 (452 _ 75^^)] | ^_ 

The formulas for ^ = 2, 3 were essentially known to Hermite. Settin g £ = 4, 5 in The- 
orem L yiel ds the less familiar and mor e complicated formulas which lEnol'skii fc Kostovl 
((199I and lEilbeck fc Enol'skiil ^9^) derived by eliminating variables from com- 
patibility conditions. Theorem L readily yields the covering map tt^ and auxiliary 
function Ke for far larger i. 
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5. Hyperelliptic reductions 



The cover tti : Ti(g2,g3) — > Eg^_g.^ introduced as part of the Hermite-Krichever 
Ansatz, i.e., the map {B^v) i— > (xojyo); is of independent interest, since explicit 
examples of coverings of elliptic curves by higher-genus algebraic curves are few, 
and the problem of determining which curves can cover Eg^ g^^ for either specified 
or arbitrary values of the invariants §2,93, remains unsolved. Ti{g2,g3) generically 
has genus g = £, as noted, and the cover will always be of degree N = £{£ + l)/2. 
The formula for xq = xo{B; 32, 53) given in Theorem L is consistent with this, since 
N equals max(A'^ + 2Nt\, Nj^ + 2NtY), the maximum of the degrees in B of the 
numerator and denominator of xq. The degrees Nj, Nj^ were computed in §|3 and 
the twisted degrees Nt\,NtY were also (see proposition I3.5|) . 

Since r^((?2,53) is hyperelliptic (defined by the irrationality ly^ ~ Le{B; g2, gs)) 
and Eg^^g^ is elliptic (defined by the irrationality = Axq — g2Xo — 53), the map tti 
enables certain hyperelliptic integrals to be reduced to elliptic ones. Just as Eg^^g^ 
is equipped with the canonical holomorphic 1-form dxo/yo, so can Ti{g2,g3) be 
equipped with the holomorphic 1-form dB/u. Any integral of a function in the 
function field of a hyperelliptic curve (here, any rational function R{B, v)) against 
its canonical 1-for m is called a hyperell iptic integral. Hyperelliptic integrals are 
classified as follows ijBelokolos et ahllflSfifl . The linear space of meromorphic 1-forms 
of the form R{B, v) dB/v, i.e., of Abelian differentials, is generated by 1-forms of 
the first, second and third kinds. These are (i) holomorphic 1-forms, with no poles; 
(ii) 1-forms with one multiple pole; and (iii) 1-forms with a pair of simple poles, 
the residues of which are opposite in sign. The indefinite integrals of (i)-(iii) are 
called hyperelliptic integrals of the first, second and third kind s. They generalize 
the three kinds of elliptic integral |Abramowitz fc StegunlllQBfil chapter 17). 

Hyperelliptic integrals of the first kind are the easiest to study, since the linear 
space of holomorphic 1-forms is finite-dimensional and is spanned by B^^^ dB / v, 
r — 1, ... ,5, where g is the genus. So there are only g independent integrals of 
the first kind. A consequence of the map tt^ : Te,{g2Tg3) Eg^ g^ is that on any 
hyperelliptic curve of the form Ti{g2,g3) there are really only g — 1 independent 
integrals of the first kind, modulo elliptic integrals (considered trivial by compar- 
ison). Changing variables in J dxa/y^, the elliptic integral of the first kind, yields 



The quantity in square brackets is rational in B, and in fact is guaranteed to be a 
polynomial in B of degree less than or equal to 17 — 1, since the left-hand integrand 
is a pulled-back version of the right-hand one, and must be a holomorphic 1-form. 
Equation H5.1|l is a linear constraint relation on the g basic hyperelliptic integrals 
of the first kind. It reduces the number of independent integrals from g to g — 1. 

The cases £ = 2, 3 of IjS.lfl may be instructive. The maps {B, v) ^ {xo,yo) were 
given in (|4.8() . H4.10|I . and the degree-(2^-|- 1) spectral polynomials L^(-; 52, .93) follow 
from table 121 li £ = 2, one obtains the hyperelliptic-to-elliptic reduction 




(5.1) 




(5.2) 
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Table 7. Polynomials specifying the holomorphic 1-forms pulled back from Eg^^g^ 



Pe{B; 92,93) 



1 1 

2 B 



3 - f92 

4 - Slg,B + i^93 

5 B^-Sfg,B^ + mig,B+Sf9i 

6 B' - Sfg,B'i + 12^g,B^ + i^glB - ^^g293 

7 r6 973 „ R4_LinS1Q„ p3 , 681373 „2 p2 2145143 „ „ p , 54071875 „2 5417685 „3 

7 « —920 ^ mViligzB H Y6~92B 5 — 5253-0 H — 53 le — 52 

Q d7 1953 „ d5 I 0001 fi„ d4 I 3335445 „2 o3 34152435 „ „ n" 



^32-0 + iU»id53£! + ^^^Tt^52-« - 2 5253-0 + ^ 

^ (_ 122490225^3 ^ 937038375 ^2 ) ^ 179425125g|g3 



where the change of variables is performed by (|4.8|l . If ^ = 3, one obtains 



where the full spectral polynomial L3(_B; 92,33) is 

n /- 06 63 „ r4 I 297„ r3 , 4185 „2 p)2 18225 „ „ r 3375 „3 , 91125 „2\ 



(5.3) 



and the cha nge of variables is performed bv M.IOII. T hese reductions were known 
to Hermite l)Konigsbergeilll878t iBelokolos et al.l 198fi() . More recently, the reduc- 



tions induced by the £ = 4, 5 coverings were worked out ijEnol'skii & Kostov 19!^ 
lEilbeck fc Enorskiilll99^ . But the reductions with £ > 5 proved too difficult to 
compute. Theorem L makes possible the computation of many such higher reduc- 
tions. 

The following proposition specifies the normalization of the pulled-back 1-form. 
It follows from the known leading-order asymptotic behavior of XQ,yQ/h' as B ^ 00. 

Proposition 5.1. For all integer £ > 1, the polynomial function Pe{B; g2, gs) '■= 
[{yo/h')^^dxQ/dB](B;g2,g3) in the hyperelliptic-to- elliptic reduction formula 

Pi{B;g2,g3)'iB _ f cIxq 



\JU,{B]g2,gz) J \J^xl - g2XQ - g^ 

where the change of variables xq = xo{B; g2, ga) is given by Theorem L, equals 
£{£+l)/4: times a polynomial Pi{B; g2, gs) which is monic and of degree£—\ inB. 

The polynomials Pi are listed in Table [3 PatP^ agree with those found by 
Enol'skii et al., if allowance is made for a difference in normalization conventions. 

A complete analysis of Lame-derived elliptic covers will need to consider excep- 
tional cases of several kinds. The covering curve rf(f;2,5'3) generically has genus 
g — £, but if the Klein invariant J — .92/(52 ~ 27(7|) is a root of one of the two Cohn 
polynomials of table 01 the genus will be reduced to £ — 1. According to conjec- 
ture l3.3l this will happen, for instance, if ^ = 2 (mod 3) and 92 = (i-e., J = 0), so 
that the base curve Eg^^g^ is equianharmonic. When g is reduced to £ — 1 in this way, 
the linear space of holomorphic 1-forms will be spanned by B^~^{B — Bq) dB/i/, 
r = — 1, where Bq is the degenerate root of the spectral polynomial; but 

l)5.1|l will still provide a linear constraint on the associated hyperelliptic integrals. 



Article submitted to Royal Society 



Lame polynomials, hyperelliptic reductions and Lame band structure 29 

Another sort of degeneracy takes place when the modular discriminant A :— 
52 — 27(;| equals zero, i.e., when J = oo. In this case Eg^^g^ will degenerate to a 
rational curve, due to two or more of ei, 62, 63 being coincident. The Lame-derived 
reduction formulas, such as H5.2I) ~ (|5.3I) . continue to apply. (They are valid though 
trivial even in the case Ci — e2 = 63, in which 32 = 93 — 0-) So these formulas 
include as special cases certain hyperelliptic-to- rafionaZ reductions. 

Subtle degeneracies of the covering map tt^ can occur, even in the generic case 
when has genus t and Eg^^g^ has genus 1. The branching structure of wi is 
determined by the polynomial of table[71 which is proportional to dxp/di?. If Pi 
has distinct roots {P^^'j^Z^ , then ne will normally have 2£ — 2 simple critical points 
on Te, of the form ±i^(*))}^~;J. However, if any i?^*) is located at a band edge, 

i.e., at a branch point of the hyperelliptic {B, j/)-curve, then , 0) will be a double 
critical point. This appears to happen when £ = (mod 3) and the base curve Eg^^g^ 
is equianharmonic; the double critical point being located at {B^v) = (0,0). Even 
if no root of Pi is located at a band edge, it is possible for it to have a double root, 
in which case each of a pair of points (i?^'^ ±1^^'^) will be a double critical point. 
By examination, this happens when £ = 4 and J = — 2^5'*/3^53. 

A few hyperelliptic-to-elliptic reductions, similar to the quadratic {N — 2) re- 
duction of Legendre and Jacobi, can be found in handbooks of elliptic integrals 
ijBvrd fc Fricd maii 195 4, §§575 and 576). The Lame-derived reductions, indexed 
by £, should certainly be included in any future handbook. It is natural to wonder 
whether they can be generalized in some straightforward way. The problem of find- 
ing the genus-2 covers of an elliptic curve was intensively studied in the nineteenth 
century, by Weierstrass and Poincare among many others, and one may reason by 
analogy with results on i — 2. One expects that for all £ > 2 and for arbitrarily 
large N , a generic Eg^^g^ can be covered by some genus-^ curve via a covering map of 
degree N. Each Lame-derived covering tt^ : r^(52,53) Eg^^g^ has N — £{£ + l)/2 
and may be only a low-lying member of an infinite sequence of coverings. Gener- 
alizing the Lame-derived coverings may be possible even if one confines oneself to 
N = £{£+ l)/2. One can of course pre-compose with an automorphism of (52, 53) 
and post-compose with an automorphism of £"^2,33 (^ modular transformation). But 
when £ = 2 a quite di fferent covering map with the same degree is known to exist 
i Belokolos et al.ll98(il) . 7r2 has two simple critical points on r2, but the other degree- 
3 covering map has a single double critical point on its analogue of Fg. Both can 
be ge neralized to include a free parameter (Burnside 1892; Belokolos & Enol'ski] 
It seems possible that when £ > 2, similar alternatives to the Lame-derived 
coverings may exist, with degree £(£ + l)/2 but different branching structures. 

6. Dispersion relations 

It is now possible to introduce dispersion relations, and determine the way in 
which the Hermite-Krichever Ansatz reduces higher-^ to £ = 1 dispersion rela- 
tions. The starting point is the fundamental multi- valued function $ introduced 
in §121 As noted, if the parametrization point (a;o,?/o) on the punctured elliptic 
curve £"^2,33 \ {(^j '^)} is over xq = B G <C, then <&(•, •; xq, yo) will be a solution of 
the £ = 1 case of the Lame equation (|2.9(l . Eg^,g^ is defined by y^ — ix^ — g2X — g^, 
so the hypothesis here is that [xq, yo) should equal {B, ±-\/4i?'^ — g2B — 53). 
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In the Jacobi form (with independent variable a), resp. the Weierstrassian form 
(with independent variable u), the crystal momentum k characterizes the behavior 
of a solution of the Lame equation under a a + 2K, resp. u ^ u + 2w. Both shifts 
correspond to motion around Eg^^g^, along a fundamental loop that passes between 
(a;, y) = (ei, 0) and (oo, oo), and cannot be shrunk to a point. (If ei, 62, 63 are defined 
by (|2.1|) . this will be because y is positive on one-half of the loop, and negative on 
the other.) By definition, the solution will be multiplied by ^ := exp[ik{2K)] = 
exp[ifc(2tt>)]. It follows from the definition (|2.11(l of $ that when i = 1, 



That is, when £ — 1 the crystal momentum is given by a complete elliptic integral. 
In the context of finite-band integration theory, this is a special case of (|4.4|l . 

It was pointed out in §0]that the spectral curve Fi that parametrizes = 1 so- 
lutions can be identified with £"^2,33 itself, via the identification {B, v) = [xq, yo/2). 
This suggests a subtle but important reinterpretation of k. In %^it was introduced 
as a function of the energy parameter, here B, which is determined only up to integer 
multiples of tt/K = tt/w, and which is also undetermined as to sign. If the presence 
of yo — 2v in Hfi.l|l is taken into account, it is clear that the £ = 1 crystal momen- 
tum, called fci henceforth, should be regarded as a function not on \ {00} 9 B, 
but rather on Fi \ {(00,00)} 9 {B^v). In this interpretation, the indeterminacy of 
sign disappears. The additive indeterminacy, on account of which ki is an elliptic 
function of the second kind, remains but can be viewed as an artifact: it is due to 
ki oc log^, where ^ is the Floquet multiplier. The behavior of fci near the puncture 
{B, 1/) = (00, 00) is easily determined. It follows from Ht).l|l that as {B, ly) (00, 00), 
i.e., {xq, yo) (00, 00), each branch of fci is asymptotic to —iv/B to leading order. 
Since B = xqtI' = 2yo have double and triple poles there, respectively, it follows 
that each branch of ki has a simple pole at the puncture. 

The multiplier ^ is a true single-valued function on the punctured spectral curve 
Fi \ {(00,00)}, and moreover is entire. One can write ^ : Fi \ {(00,00)} ^ \ 
{0,00}, since the multiplier is never zero. Like ki, this function is not algebraic: 
it necessarily has an essential singularity at the puncture. The {{B, v), ^)-curve over 
Fi \ {(00,00)} 3 {B,v), which is a single cover, and the (i3,^)-curve over C 3 B, 
which is a double cover, are both transcendental curves. 

The crystal momentum for each integer i > 1 may similarly be viewed as an 
additively multi- valued function on the punctured spectral curve F^ \ {(00,00)}. 
It will be written as ki{B, g2, gs), with the understanding that for this to be well- 
defined, B, V must be related by the spectral curve relation = Li{B; g2, gs). The 
quantity ki{B,v; g2, ga) will not be undetermined as to sign. Suppose now that the 
projections tt^ : F^ ^ Eg^_g^ of the Hermite-Krichever Ansatz are regarded as maps 
ne, : F| ^ Fi. That is, -Kg maps [B, u) S F^ to the point {B' , v') := {xo,yo/2) S Fi. 
The reductions {B, v) ^ {B', v') for ^ = 2, 3, for example, follow from (|4.8|) - (|4.11() . 

Proposition 6.1. If the integration of the Lame equation on the elliptic curve Eg^^g^. 
for integer £ > 1, can be accomplished in the framework of the Hermite-Krichever 
Ansatz by maps Hi : Ti ^ Fi and K£ : F^ — > P^, where tti and ki map the point 
{B,v) to {Bi{B;g2,g3),iye{B,iy;g2,g3)) and ki{B; g2, gs)^, respectively, then the 
dispersion relation for the Hermite-Krichever solutions will be k — ki{B, i^; (72,33)7 
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in which can he expressed in terms of ki by 

ke{B,i^;g2,g3) = ki{Bt{B,v;g2,gs.),vt{B,v\g2,g3)) -ike,{B,v;g2,gz)v. (6.2) 

This proposition follows immediately from the form of the Hermite-Krichever 
solutions (|4.1|l . The first term in Ht).2|) arises from the <i>, $' factors, and the second 
from the exponential. The factors A, B in H4.1|l do not contribute to the crystal 
momentum. Equation Ht).2|) could also be derived from the general theory of finite- 
band integration, specifically from the formula (|4.4() . However, a derivation from 
the Hermite-Krichever Ansatz seems more natural in the present context. 

The effort expended in replacing two-valuedness by single- valuedness is justified 
by the following observation. As a function of B alone, rather than of the pair (_B, v), 
each of the two terms of (|6.2|) would be undetermined as to sign. This ambiguity 
could cause confusion and errors. The present formulation, though a bit pedantic, 
facilitates the determination of the correct relative sign. 

The form of the (S, v) ^ [B' , v') map assumed in the proposition is of course 
the form supplied by Theorem L. Substituting the formulas of the theorem into (|6.2(l 
yields an explicit expression for ki. It follows readily from this expression that for all 
integer ^ > 1, each branch of k^ on the hyperelliptic (-B, i/)-curve Ti{g2, gs) satisfies 
ke ~ —w/B^, {B, v) — > (oo, oo). Since B. v have double and order-(2£+ 1) poles at 
the puncture (oo,oo), respectively, this implies that irrespective of each branch 
of the crystal momentum has a simple pole at the puncture. 

7. Band structure of the Jacobi form 

The results of §§13 through IHl were framed in terms of the elliptic-curve algebraic 
form of the Lame equation. Most work on Lame dispersion relations has used the 
Jacobi form instead, and has led accordingly to expressions involving Jacobi theta 
functions. To derive dispersion relations that can be compared with previous work, 
the formulation of §|nimust be converted to the language of the Jacobi form. 

The relationships among the several forms were sketched in §13 In the Weier- 
strassian and Jacobi forms, the Lame equation is an equation on C with doubly 
periodic coefficients, rather than an equation on the curve £"92,33- I'^ the conver- 
sion to the Jacobi form the invariants 52 , 53 are expressed in terms of the modular 
parameter m by H2.2|l . with A = 1 by convention. The coordinate x on Eg^^g^ 
is interpreted as the function p{u; g2, ga), i.e., as msn^{a\m) — i(m + 1), where 
u G C and a := u — iK'(TO) € C are the respective independent variables of 
the Weierstrassian and Jacobi forms. The holomorphic differential dx/y corre- 
sponds to du or da, and the derivative yd/dx to d/du or d/da. The coordi- 
nate y = {yd/dx)x on Eg^^g^ is interpreted as p'(u; (72, ffs), i-e., as the doubly peri- 
odic function 2m sn(a|m) cn(a|m) dn(a|m) on the complex a-plane. The functions 
y/x — e-y, 7 = 1, 2, 3, correspond to — idn(Q:|rn), — i m^/^cn(Q;|m) and m^^'^ sn{a\m). 

B = -E+l£{£+l){m+l) (7.1) 

relates the accessory parameters B,E of the different forms. 

With these formulas, it is easy to convert the Lame polynomials of table [5] to 
polynomials in sn(a|TO), cn(a|TO), dn(Q;|m), and the spectral polynomials^ of table|3l 
to polynomials in E, for comparison with the list given bv Rrscottl ljl964l §9.3.2). 
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Definition. The Jacobi-form spectral polynomial L^(£'|to) is the negative of the 
spectral polynomial Li{B; g2, gs), when B;g2,g3 are expressed in terms of E,m. 
It is a monic degree- (2i? + 1) polynomial in E with coefficients polynomial in m, 
and can be regarded as Jlslol-^ ^ Es{m)], where the roots {Eg} are the values 
of the energy E for which there exists a Lame polynomial solution of the Lame 
equation, counted with multiplicity. (The negation is due to the relative minus sign 
in the B ^ E correspondence l|7.1fl .) 

Definition. The £'th Jacobi-form spectral curve :— ^({171) is the hyperelliptic 
curve over 9 i5 comprising all {E,i>) satisfying — Le{E\m). In the usual 
case when m is real, i' will be real if _B is in a band, and non-real if is in a 
lacuna. In both cases it is determined only up to negation. By convention, the 
correspondence between the curve 9 (E, v) and the previously introduced curve 
F^ 3 {B, v), which was defined by v'^ = \-i{B\ g2, 53), is given hy v = +0). 

The following cases are examples. When £ = 1,2,3, the spectral polynomial 
factors over the integers into polynomials at most quadratic in E. In full, 

Li{E\m) ^ {E - 1){E - m){E - m - 1) (7.2) 

l2{E\m) ^ [E^ - 4(m + 1)E + 12m] [E - m - 1){E ~ 4m - 1){E ~ m ~ 4) 

(7.3) 

L3{E\m) = {E-4m-4) [E^ - 2(2™ + 5)E + 3(8m + 3)] (7.4) 
X [E^ - 2(5m + 2)E + 3(3™^ + 8m)] 
X [E^ - 10(m + 1)E + 3(3™^ + 26m + 3)] . 

In (|7.3|) and (|7.4|) the first factor arises from L\{B; g2, ga) and the remaining three 
from the factors L^iB; e^, (72, 53), 7 = 1,2, 3. In (|7.2I) there is no Type I factor. The 
polynomials (|7.2I) - (|7.4I) agree with those obtained by Arscott. 

The derivation of the Jacobi-form s pectral polyn omial Li{E\m) from (i?; g2, 53) 
sheds light on a regularity noticed by lined |)1940qI § 7) , which arises in the lemnis- 
catic case m — ^. Ince observed that if £ < 6, at least, then L^(£'|i) has an integer 
root, namely E — £{£ + l)/2. In fact, this is the case for all integer £. By 1)7. 
the presence of this root is equivalent to the full spectral polynomial L£{B]g2,0) 
having i3 = as a root. But if m = i, it follows from H2.1|l that 62 = 0. A glance at 
the pattern of coefhcients in table |2| reveals that if 93 = and a singular point e-y 
also equals zero, then either the Type-I spectral polynomial L\{B; g2, 53) (if £= 0,3 
mod 4) or one of the three Type-II spectral polynomials LfiB; e^, 32, 53) (if £= 1,2 
(mod 4)) will necessarily have B = as a root. 

Dispersion relations in their Jacobi form can now be investigated. Recall that if 
£ = 1, the Jacobi-form Lame equation Hl.l() has $(•; aol"^') as a solution, where the 
theta quotient $ (the Jacobi-form version of Halphen's relement simple) is defined 
in (|1.3|) . and the multi- valued parameter ao is defined by dn^(ao|"^) = E — m. 
This solution has crystal momentum k = ki equal to — iZ(ao|m)-|-7r/2K(m), which 
is undetermined up to sign, and is also determined only up to integer multiples 
of 7:/K{m). The sign indeterminacy is due to dn^(-|m) being even. This causes ao 
to be undetermined up to sign, and ki as well, since the function Z(-|m) is odd. 

The parametrization point ao, or the equivalent point uq := ap -I- iK'(m) of the 
Weierstrassian form, corresponds to the parametrization point {xq, ?/o) of the funda- 



Article submitted to Royal Society 



Lame polynomials, hyperelliptic reductions and Lame band structure 33 



mental function on the elliptic curve Eg^^g^. The correspondence is the usual one: 
xq = p{uq] (72, ffs), yo = p'{uo, 92,93)- The first of these two equations says that xq — 
m sn^(Q!o|77i) — ■|('Ti + 1), and the latter that = 2m sn(Q!o|m) cn(Q;o|m) dn{aQ\m). 

The formula which computes ao from E, namely dii' {ao\m) = E — m, is readily 
seen to be a translation to the Jacobi form of the familiar condition a;o — B, which 
simply says that the parametrization point (a;o,?/o) must be 'over' B € C. 

The correspondence between the Jacobi and elliptic-curve forms motivates the 
following reinterpretation of the crystal momentum of the fundamental solution $, 
which is modelled on the reinterpretation of the last section, ki should be viewed as 
a function not of the energy E € C, but rather of a point {E, v) on the punctured 
Jacobi-form spectral curve Ti{m) \ {(00,00)}. There are two such points for each 
energy E, except when £^ is a band edge. This is the source of the sign ambiguity 
in the parameter ao. Since yo = 2v = 2ij>, the equation 

m sn(ao|m) cn(ao|TO) dn(Q!o|m) — ii> (7.5) 

determines a unique sign for ao, provided that v is specified in addition to E. ki will 
be written as ki{E, v\m), with the understanding that for this to be well-defined, 
the pair E,v must be related by the spectral curve relation iP — Li{E\ra). The 
additively undetermined quantity ki{E,v\m) will not be undetermined as to sign. 
It is easily checked that on each branch, ki ^ v/ E as {E, v) — (00, 00). 

Definition. A solution of the Jacobi-form Lame equation is said to be an 

Hermite-Krichever solution if it is of the form 

^(sn^(a|m)) $(a; ao|m) -f 2S (sn^(a|TO)) $'(a; aolm) exp(Ka), (7.6) 

for some ao G C and k G C. Here A, B are polynomials, and $' := (d/da)$. 

The expression (|7.6|l is a replacement for the original Jacobi-form expression 
()1.5(l . to which it is equivalent. Regardless of which is used, it is easy to compute 
the crystal momentum of an Hermite-Krichever solution. The momentum computed 
from 1)7.6(1 will be [— i Z(ao|m)-|-7r/2K(m)] — Ik, up to additive multi-valuedness. The 
first terrn arises from the factors, and the second from the exponential. The 

factors A,B do not contribute, since sn^(a|m) is periodic in a with period 2K(m). 

The Jacobi form of the Hermite-Krichever Ansatz asserts that for all integer £ 
and TO e C \ {0, 1}, there is an Hermite-Krichever solution for all but a finite 
number of values of the energy E. On the elliptic curve Eg^ ^g^ , these solutions were 
constructed from two maps: a projection tt^ : — > Eg^_g.^ and an auxiliary function 

: ^ P^. But TTg should really be regarded as a map from to Fi, on account 
of the correspondence between Eg^^g^ and Fi provided by (xo,?;o) = {B,2v). The 
following is the Jacobi-form version of proposition l6.1l 

Proposition 7.1. Suppose that the integration of the Lame equation on the ellip- 
tic curve Eg^^g^, for integer i > I, can be accomplished in the framework of the 
Hermite-Krichever Ansatz by the maps ni : Ti ^ Fi and ki : F^ — > P^, where 
and Kt map the point {B,v) to (5^(5; 32, 53), t'f (S, 32, ffa)) and ki{B] g2, 93)^ , 
respectively. Then the dispersion relation for the solutions of the Jacobi form of the 
Lame equation will be k ~ ki(E,i>\m), where up to additive multi-valuedness 

ki{E, i?\m) = fci (^£e{E\m), i>i{E, i>\m) | m) -I- ki{E\m)i>, (7.7) 
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k/(V2K) 

Figure 1. Dispersion relations for £ = 1, 2, 3 in the lemniscatic case m — ^. 



in which 



£i{E\m) 
vt^E, i/|m) 
kt{E\m) 



= -Bt[-E 
= ki{-E+\l{l 



+ \){m + l);52(TO),ff3("i)) + |(m 
+ l)(m + 1), g2(jn), gsim)) , 
l){m + 1): g2{m), gsim)) . 



1), 



(7.8) 
(7.9) 
(7.10) 



The formula H7.7|l follows by inspection. The projection ni reduces the integra- 
tion of the Lame equation to the integration of an £ = 1 equation, the 'B' parameter 
of which equals Be{B; g2, 53). By 17.111 . the '£" parameter of the 1=1 equation will 
be the right-hand side of (|7.8I) . The two terms of H7.7|l are simply the two terms of 
[— iZ(Q;o|'7i) -I- 7r/2K(m)] — Ik. The equality = it> has been used. 

It is straightforward to apply Droposition l7. ll to the cases £ = 2,3, since the cov- 
erings 7r2 , TTa and auxiliary functions K2 , K3 were worked out in §0] A brief discussion 
of the £ = 2 case should suffice. After some algebra, one finds 



£2{E\m) 



E^ - 12(m + ifE - 4(to + l){Am'^ - 19m + 4) 
9[E^ - 4(m + 1)£; + 12m] 



, , (£; + 2m-4)(£;-4m-|-2)(i;-4m-4) 
iy2(E, v\ni) = — u. 



K2(£'|m) 



27 [£:2 - 4(m + 1)£; + 12m]' 
2 



2,[E^ - ^{711 + l)E + I2m\ 



(7.11) 
(7.12) 
(7.13) 



from which k2{E,i>\m) may be computed by (|7.7|l . Like ki, k2 is determined only 
up to integer multiples of tt/K := n/K{m). Each branch of has the property 
that k2{E, i/\m) ~ ±£'^/^, E +00, with '±' determined by the sign oi v = v{E). 
This is a special case of a general fact: for all integer 1>1, ki{E,i'\m) ^ ±i?^/^, 
E — > -l-oo, since each branch of ki is asymptotic to {—Y^^v/E^ as {E, v) {00, 00). 

The real portions of the dispersion relations k = ki(E, i>\^), k2{E,i>\^) and 
k^[E, i/|i) are graphed in figure ^ For ease of viewing, each crystal momentum is 
regarded as lying in the interval [0, i: /2K]; which is equivalent to choosing the sign 
of i> = ^{E) in an iJ-dependent way. As 17.2|l - 17.4|l imply, the two 1=1 bands are 

1], [|, 00), the three £ = 2 bands are [3 ^ VS, |], [3, |], [3 -I- -s/S, 00), and the four 
£ = 3 bands are [| - ^6,6- Vl5], - V6,6], [| + v%,6-(- VTS], -f V6,oo). 
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The £ = 1 graph ag rees with that of i Li et al.1 l)2000l figure 6), and for con- 
firmation, with that of ISutherlandl l)l973l figure 1). Unfortunately, the £ — 2 
graph disagrees with that of Li et al. in the placement or direction of curvature 
of each of the two upper bands. The algorithm they used for reducing £ — 2 
to £ = 1, which was based on Hermite's solution of the Jacobi-form Lame equation 
l|Whittaker fc Watsonlll927l §23.71), evidently yielded incorrect results for these 
bands. It appears that for the middle band, at least, the discrepancy can be traced 
to an incorrect choice of relative sign for the two terms of fc = A:2- The reinterpre- 
tation of the crystal momentum as a function on the spectral curve, rather than a 
function of the energy, eliminates such sign ambiguities. 



8. Summary and final remarks 

A new approach to the closed-form solution of the Lame equation has been in- 
troduced. Theorem L provides a formula for the covering map of the Hermite- 
Krichever Ansatz in terms of certain polynomials which are of independent inter- 
est, namely twisted spectral polynomials. The theorem permits an efficient com- 
putation of Lame dispersion relations, of a mixed symbolic-numerical kind. Cohn 
polynomials, which are a new concept, have also been introduced. The roots of 
such a polynomial are the points in elliptic moduli space at which a Lame spectral 
polynomial has a double root, so that the Lame spectral curve becomes singular. 
Twisted and theta-twisted Cohn polynomials could be defined, as well. 

The approach of this paper can be extended from the Lame equation to the 
Heun equation, which as a differential equation on the elliptic curve Eg^^g^ has up to 
four regular singular points, positioned at the finite Weierstrass points {(e^, 0)}^^^ 
as well as at (oo , oo). Its Weierstrassian form is called the Treibich-Verdier equation 
llbmirno vl2002l) . and its Jacobi form, at least when only two of the Weierstrass points 
are singular points, the associated Lame equation ( Magnus fc Winkle^llQTfll § 7.3). 



The 'four triangular nu mbers' condition for the H eun e quation to have the 
finite-b and property, due to iTreibich fc Verdieil l)l992l) and lOesztesv fc WeikardI 
I l995fil) . IS now well known. In the finite-band case, the nu mber of points in the 



algebr aic spectrum has been computed up to multiplicity ijGesztesv fc WeikardI 
Il995fefl . The corresponding band-edge solutions are Heun polynomials. Applying 
the Hermite-Krichever Ansatz to the finite-band Heun equation leads to a greater 
variety of coverings of -Ega.sa than arise in the solution of the integer-f Lame equa- 
tion; for example, coverings b y a genus-2 hyperelliptic curve that have degrees 3, 4, 5 
l)Belokolos fc Enorskiilf2000|) . These coverings play a role in the construction of el- 
liptic s oliton solutions of certai n nonlinear evolution equations that occur in fibre 
optics ( Christiansen et aT]l20nr^ . A treatment of the Heun equation along the lines 



of this paper will appear elsewhere. 
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